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ABSTRACT Microarrays containing 1046 human cDNAs 
of unknown sequence were printed on glass with high-speed 
robotics. These 1.0-cm 2 DNA "chips" were used to quantita- 
tively monitor differential expression of the cognate human 
genes using a highly sensitive two-color hybridization assay. 
Array elements that displayed differential expression patterns 
under given experimental conditions were characterized by 
sequencing. The identification of known and novel beat shock 
and phorbol ester-regulated genes in human T cells demon- 
strates the sensitivity of the assay. Parallel gene analysis with 
microarrays provides a rapid and efficient method for large- 
scale human gene discovery. 



Biology has entered the genome era (1). Complete genome 
sequences for all of the model organisms and human will 
probably be available by the year 2003 (2). Torrents of human 
expressed sequence tags (ESTs) provide a starting point for 
elucidating the function of tens of thousands of cognate genes 
(3). Genome analysis will provide insights into growth, devel- 
opment, differentiation, homeostasis, aging, and the onset of 
diseases ( 1-3 ). A detailed understanding of the human genome 
will require the implementation of sophisticated methods for 
gene expression analysis and gene discovery. 

Recently, a microarray-based method for high-throughput 
monitoring of plant gene expression was described {A). This 
"chip"-based approach involved using microarrays of cDNA 
clones as gene-specific hybridization targets to quantitatively 
measure expression of the corresponding plant genes (4, 5). A 
two-color fluorescence labeling and detection scheme facili- 
tated sensitive differential expression analysis of different 
plant tissues (4, 5). The efficiency of this approach for studies 
in higher plants suggested the use of this method for human 
genome analysis (4-7). Here, we report the use of cDNA 
microarrays for human gene expression monitoring, biological 
investigation, and gene discovery. 

MATERIALS AND METHODS 

Human cDNA Clones. The cDNA library was made with 
mRNA from human peripheral blood Ivmphocvtes trans- 
formed with the Epstein-Barr virus. Inserts >600 bp were 
cloned into the lambda vector AVES-R to generate 10 7 -10" 
recombinants. Bacterial transformants were obtained by in- 
fecting £. colt strain JM107/AKC. Colonies were picked at 
random and propagated in a 96-well format, and minilvsate 
DNA was prepared by alkaline lysis using REAL preps 
(Qiagen, Chatsworth. CA). Inserts were amplified bv PCR in 
a 96-well format using primers (PAN132, 5 : -CCTO 
TATACTTTA ACGTCA AGG ; and PAN133, 5'-TTGTGTG- 
GAATTGTGAGQGG) complementary to the AYES 
polylinker and containing a six-carbon amino modification 
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(Glen Research, Sterling. VA) on the 5* end. PCR products 
were purified in a 96-well format using QlAquick columns 
(Qiagen). 

Microarray Preparation. Amino-modified PCR products 
were suspended at a concentration of 0.5 mg/ml in 3x 
standard saline citrate (SSC) and arrayed from 96-well micro- 
liter plates onto silylated microscope slides (CEL Associates. 
Houston) using high-speed robotics (4-7). A total of 1056 
cDNAs. representing 1046 human clones and 10 Arabidopsis 
controls, were arrayed in 1.0-cm : areas. Printed arravs were 
incubated for 4 hr in a humid chamber to allow rehydration of 
the array elements and rinsed, once in 0.2% SDS for 1 min, 
twice in H;0 for 1 min, and once for 5 min in sodium 
borohydride solution (1.0 g of NaBH4 dissolved in 300 ml of 
PBS and 100 ml of 1009c ethanol). The arravs were submerged 
in H : 0 for 2 min at 95°C, transferred quickly into 0.2% SDS 
for 1 min. rinsed twice in H : 0, air dried, and stored in the dark 
at 25°C. 

Fluorescent Probes. Tissue mRNAs were purchased 
(CLONTECH). Jurkat mRNA was isolated as described by 
Schena et aL (4). Probes were made as described (4) with 
several modifications. The reverse transcriptase used here was 
Superscript II RNase H- (GIBCO). The Cy5-dCTP was 
purchased from Amersham. Each reverse transcription reac- 
tion contained 3.0 *ig of total human mRNA." Arabidopsis 
control mRNAs were made by in vitro transcription of cloned 
H AT4, HAT21 and YesAt-23 cDNAs (4. 8. 9) using an RNA 
Transcription Kit (Stratagene). For quantitation, the mRNAs 
were doped into the reverse transcription reaction at ratios of 
1:100.000. 1:10.000. and 1:1000 (wt/wt) respectively. Following 
the reverse transcription step, samples were treated with 2.5 u\ 
of 1 M sodium hydroxide for 10 min at 37°C. then neutralized 
by adding 2.5 u\ of 1 M Tris HCI (pH 6.8) and 2.0 M l of 1 M 
HC1. Probe mixtures contained cDNA products derived from 
3 jig of total mRNA. suspended in 5.0 u.1 of hybridization 
buffer (5x SSC plus 0.2% SDS). 

Hybridization and Scanning. Probes were hybridized to 
1.0-cm : microarrays under a 14 x 14 mm glass coverslip for 
6-12 hr at 60°C in a custom-built hybridization chamber (4-7). 
Arrays were washed for 5 min at room temperature (25°C) in 
low stringency wash buffer (lx SSC/0.2% SDS), then for 10 
min at room temperature in high stringency wash buffer (0.1 x 
SSC/0.2% SDS). Arrays were scanned in 0.1 x SSC using a 
f luorescence laser scanning device (4-7). fitted with a custom 
filter set (Chroma Technology. Brattleboro, VT). Accurate 
differential expression measurements (i.e., final fluorescence 
ratios ) were obtained by taking the average of the ratios of two 
independent hybridizations. 



Abbreviation: EST, expressed sequence lag. 

Data deposition: The sequences reported in this paper have been 
deposited in the GenBank database (accession nos. U56654-U56660). 
n*o whom reprint requests should be addressed, e-mail: schenafo 
cmgm.stanford.edu. 
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Cell Culture. Jurkai cells were grown in a tissue culture 
incubator (3TC and 5Sc CO:) in RPM1 medium supplemented 
with 10Tc fetal bovine serum. 100 u.g of streptomycin per ml. 
and 500 units of penicillin per ml. Heat shock corresponded to 
a 4-hr incubation at 43 C C. Phorbol ester treated cells were 
grown for 4 hr in the presence of 50 ng of phorbol 12-myristate 
13-acetate (PMA) per ml. 

RNA Blotting. Dot blots were performed as described (4). 

DNA Sequencing. Sequences were obtained using the 
PAN132 and PAN133 primers and a 373A automated se- 
quencer, according to the instructions of the manufacturer 
(Applied Biosystems). 

Computer Graphics and Informatics. Pseudocolor represen- 
tations of fluorescent images were made with National Institutes 
of Health image software (version 1.52). Software for differential 
expression representations was purchased from Imaging Re- 
search (St. Catherine s. ON. Canada). Sequence searches were 
made to the nonredundam nucleotide data base at the National 
Center for Biotechnology Information (NCBI) using Macintosh 
blast software. The EST data base was accessed via the World 
Wide Web (hup: /www.ncbLnlm.nih.gov/ ). 

RESULTS 

Gene Discovery and the Heat Shock Response. MicToarrays 
were used to examine the heat shock response in cultured 
human T (Jurkai) cells. Control (37°C) and heat-treated 
(43°C) cells were harvested and lysed. and total mRNA from 
the two cell samples was labeled by reverse transcriptase 
incorporation of fluorescein- and Cy5-dCTP. respectively. In 
a second set of labeling reactions, the fluorescent groups were 
"swapped" such that samples from control and heat-treated 



samples were labeled with Cy5- and f luorescein-dCTP. respec- 
tively. Each pair of fluorescent probes was hybridized to a 
1056-element microarray. The arrays were washed at high 
stringency and scanned with a confocal laser scanning device 
to detect emission of the two fluorescent croups. 

Hybridization signals were observed to >95r> of the human 
cDNA array elements, but not to any of the Arabidopsis 
negative controls (Fig. 1). Fluorescence intensities spanned 
more than three orders of magnitude for the 1046 array 
elements surveyed (Fig. 1 ). Comparative expression analysis of 
heat shocked versus control cells in the two experiments 
revealed 17 array elements that displayed altered fluorescence 
ratios of £2.0-*fold (Figs. 1 and 24). Of the 17 putative 
differentially expressed genes. 11 were induced by heat shock 
treatment and 6 displayed modest repression (Figs. 1 and 2.4). 

To determine the identity of the heat-regulated genes. 
cDNAs corresponding to each of the 17 array elements were 
sequenced on the proximal and distal end. Data base searches 
revealed perfect matches for 14 of the 17 clones, and in each 
case proximal and distal cDNA sequences mapped to the same 
gene (Table 1). Of the 1046 human genes examined on the 
microarray. the five most highly induced in heat-treated cells 
were heat shock protein 90a (hsp90o), dnaJ, hsp9O0. poiyu- 
biquitin. and t-complex polypeptide- 1 (tcp-1 ) (Table 1). Three 
of the 17 clones did not match any entry in the public data base, 
though one of the clones (B7) exhibited significant homology 
to an EST from Caenorhabditis elegans (Table 1). Each of the 
novel sequences (B7-B9) exhibited ~2-foId induction (Table 1) 
and relatively low-level expression (Table 2). 

To confirm the microarray results. mRNA levels for each of 
the genes were measured by RNA blotting. Each of the genes 
that displayed heat shock induction, including the three novel 



-Hut Shock 



♦Heat Shock 



•□»»■• * 



. » * . rn» -» - o*> o ; ■» , 

-0^» 09 -> m 

OQ o * c 

3 © O 33 *30> O 

333 [>[L O --0#33 33 

3 - >^\SP 9 '3) 3 O 33 

3 39 • 3 ^ i > _ O 3 

- -o& M o o - -> j 

O 3 , - O 

» * - • Cf3]o 3 » 3 ><* * » 

»: • j - * m » - * .. • • * - 

* ^ ^ „ * + > ^ » - 
*• * 9m 9 » »[•] » — 

a«^3. _ » i * -> > m ■ 

* • - j i . * . 
fy*] » *..»»•• j » 

• # * • « ft 9 * - • * • * 

• ■ r~i • * ■ • * * • • 




10 100 
Expression Level (per 100,000) 



Fig. 1. Human gene expression monitored on a microarray. Fluorescent scans represented in a pseudocolor scale correspond to expression leveb. 
The array contains \Q Arabidopsis controls (upper left corner, elements 1-10) and 1046 human peripheral blood cDNAs Fluorescent probes were 
prepared by labeling mRNA from Jurkat cells grown at 37X (-Heat Shock, A) or 43°C ( +Heat Shock. B). Array elements that display altered 
fluorescence intensity (white boxes) corresponded to genes activated (red boxes) or repressed (green boxes) by heat shock. The color bar was 
calibrated in separate experiments using known quantities (wt/wi) of Arabidopsis control mRNAs added to the labeling reaction. Microarray rows 
(at left) and columns (at the top) are demarcated at 10 clement increments (white circles). (Bar « 1 mm.) 
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Expression Ratios 

Fig. 2. Elemental displays of activated and repressed genes. Fluorescence ratios of two-color microarrav scans (Fig. 1) are depicted 
schematically. Fluorescein-labcled probes from Jurkat cells subjected to (A) heat shock or (£) phorbol ester treatment wire compared with 
Cy.vlabeled probes from untreated cells. In a second set of reactions, the fluorescent groups were swapped (sec text). The data represent the average 
of the ratios from two hybridizations, excluding values in which the difference of the two ratios was greater than half the average ratio. The color 
bar corresponds to expression ratios, which are independent of the absolute expression level of a erven gene. 



Table 1. MicToanay elements corresponding to differentially expressed genes 



Clone • 


Row 


Column 


Ratio 


Blast identity 


Accession no. 


Bl 


24 


21 


0.5 


CYC oxidase III 


J01415. J01415 


B2 


1 


31 


0.5 


0-Actin 


NR, X00351 * 


B3 


15 


8 


0.5 


CYC oxidase III 


J01415. J01415 


B4 


32 


19 


0.5 


CYC oxidase III 


J01415, J01415 


B5 


17 


8 


0.5 


CYC oxidase III 


J01415. J01415 


B6 




31 


0.5 


0-Actin 


NR. X00351 


B7' 


5 


4 


2.0 


Novel* 


U56653. U56654 


BS 


2 


19 


2.0 


Novel * 


U56655, U56656 


B9 


14 


5 


^ i , 


Novel* 


t'56657. U56658 


BIO 


7 


8 


2.4 


Polvubiquitin 


X04803. X04K03 


Bll 


i: 




2.4 


TCP-1 


X52K82. X52R82 


bi: 


28 




2.5 


Polvubiquitin 


M17597. M17597 


B13 


14 


7 


2.5 


Polvubiquitin 


X04803. X04H03 


B14 


20 


9 


2.6 . 


HSP9O0 


Ml 6660. Ml 6660 


B15 


30 


12 


4.0 


DnaJ homoloc 


D133RK. D133K8 


B16 


10 


5 


5.8 


HSP90o 


X07270. X07270 


B17, 


13 


16 


6.3 


HSP9(»o 


M27024. X15183 


B18 


7 


19 


2.0 


^-microglobulin 


S54761. M30683 


B19 


21 


30 


2.1 


Novel 1 


U56659. U56660 


b:o 


3 


26 


2 2 


/^microglobulin 


S54761. M30683 


b:i 


1 


18 


2.6 


PGK 


Ml 1968, L00160 


b:: 




30 


3.5 


NF-«B1 


Z47744, M55643 


B23 


20 


16 


19 


PAC-1 


LI 1329. L11329 



Clone name, array position (Fig. 1). fluorescence ratio, sequence identity, and accssion number of cDNAs that manifested 
a differential expression pattern with probes prepared from heat shock- (Bl-17) or phorbol ester-treated (B18-23) Jurkat cells. 
Clones showing >98<£ identity over 300 nucleotides were assumed to be identical to known sequences. All genes are nuclear 
except CYC'oxidase III (mitochondrial). Accession numbers reflect the highest score for proximal and distal sequence traces, 
respectively. CYC. cytochrome e: TCP-1. T-complex polypeptide; HSP, heat shock protein: PGK. phosphoglvccratc kinase; 
NF-kB, nuclear factor- kappaB; PAC-1. phosphatase of activated cells: and NR, trace not readable due to* the presence of 
poly(A)+ Iran. 

*B7 is 679e identical to an EST from C. eleganj (D76026). 
f No match in the public data bases. 
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Table 2. Human gene expression monitored by microajray and RNA blot analyses 
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Gone 



Expression leveL per 10-* mRXAs 



Blast identitv 



Bl CYC oxidase III 

B2 0-Actin 

B3 CYC oxidase III 

B4 CYC oxidase III 

B5 CYC oxidase III 

B6 0-Actin 

B7 Novel (wcakJv to D76026) 

B8 Novel 

B9 Novel 

BIO Polvubiquitin 

Bll TCP-1 

B12 Polvubiquitin 

B13 Polvubiquitin 

B14 HSP9O0 

B15 DnaJ homolog 

B16 HSP90o 

B17 HSP90o 

B18 fc-microglobulin 

B19 Novel 

B20 ^microglobulin 

B21 Phosphogrvcerate kinase 

B22 NF-KBl 

B23 PAC-1 



Microarray 


Rjuo 


RNA blot 


Ratio 


92/46 


03 


100 /SO 


0.8 


240/120 


05 


270/280 


1.0 


36/18 


OS 


ND 


ND 


76/38 


OS 


ND 


ND 


62/31 


OS 


ND 


ND 


180/89 


03 


ND 


ND 


13/2.6 


2.0 


0.77/1.8 


•» ^ 


2.0/4.0 


2.0 


13/3.4 


"» i 


0.8/1.8 


;_2 


1.2/1.8 


t t 


0.8/1.9 


2.4 




J.O 


23/5.5 


2.4 


7.1/27 


3.8 


0.8/2.0 


23 


ND 


ND 


1.7/43 


23 


ND 


ND 


75/200 


2.6 


30/120 


4.0 


i:o/4.o 


4.0 


1.6/13 


8.1 


0.6/33 


5.8 


3.2/29 


9.1 


0.8/5.0 


63 


8.6/62 


7.2 


1.0/2.0 


2.0 


5.4/15 


2.8 


1.2/23 


2.1 


43/93 


23 


2.7/5.9 


22 


ND 


ND 


2.4/6.2 


2.6 


4.7/9.2 


2.0 


1.7/6.0 


33 


0.65/4.7 


7.2 


0.5/9.5 


19 


0.21/15 


71 
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d t m d untreated cells. Clone and gene names are given in Table 1. ND. not 



sequences, exhibited elevated mRNA levels by dot blot analysis 
(Table 2). In all cases, expression ratios as determined by "the 
two procedures differed by <2-fold for the genes identified in 
the heat shock experiments (Table 2). The two assavs differed 
more widely in terms of assessing absolute expression levels; 
nonetheless, absolute expression as monitored on a microarray 
typically correlated with RNA blots to within a factor of five 
(Table 2). 

Phorbol Ester Signaling. To explore a signaling pathway 
distinct from the heat shock response. microatTays were used 
to examine the cellular effects of phorbol ester treatment. 
Jurkat cells were treated with phorbol ester, harvested, lysed, 
and used as a source of mRNA. Samples of mRNA from 
untreated or phorbol ester-stimulated cells were labeled with 
reverse transcriptase. The probes were mixed, hybridized to 
microarrays. and scanned for fluorescence emission of the two 
fluorescent groups. A total of six array elements displayed 
s=2.0-fold elevated signals with probes from phorbol ester- 
treated cells relative to control samples (Fig. 2B). 

To determine the identity of the phorbol ester-induced 
genes, clones corresponding to the six array elements were 
sequenced. Data base searches revealed perfect matches for 
five of the six sequences (Table 1). The two most highlv 
induced genes were the PAC-I tyrosine phosphatase and 
nuclear factor-kappa Bl {NF-kB1)\ modest activation was 
observed for phosphoglycerate kinase and ^-microglobulin 
(Table 1). One remaining clone (B19) did not match any entrv 
in the public data base (Table 1). B19 displayed a 2.1-fold 
induction and. similar to the novel heat shock genes, a rela- 
tively low absolute expression level (Tables 1 and 2). All six of 
the phorbol ester-inducible genes displayed increased steady- 
state mRNA levels tjy RNA blotting (Table 2). PAC-1 expres- 
sion (Fig. 1; Table 2) defined a detection limit of -1:500.000 
for the assay. 

Transcript Imaging in Human Tissues. To determine 
whether microarrays could be used to monitor expression in 
human tissues, probes were prepared from human bone mar- 



row, brain, prostate, and heart by labeling each mRNA sample 
with Cy5-dCTP. In a separate reaction, a control probe was 
prepared by labeling Jurkat mRNA with fluorescein-dCTP. 
The four Cy5-labeled probes were each mixed with an aliquot 
of the fluorescein-labeled control sample, and the four mix- 
tures were hybridized to separate microarrays. The arrays were 
washed and scanned for fluorescence emission, and hybrid- 
ization signals for each of the tissues samples weft normalized 
to the Jurkat control to generate an expression profile for each 
of the 1046 clones present on the array. 

Delectable expression was observed' for all 15 of the heat 
shock and phorbol ester-regulated genes in the four, tissue 
types examined (Fig. 3). In general, the expression level of each 
gene in Jurkat cells correlated rather closely with expression in 
the four tissues (Table 2: Fig. 3). Genes encoding 0-actin and 
cytochrome c oxidase, the two most highly expressed of the 15 
genes in Jurkat cells (Table 2), were highly expressed in bone 
marrow, brain, prostate, and heart (Fig. 3A). Expression of 
cytochrome c oxidase. hsp90o. and the novel B7 sequence was 
significantly greater in heart than in the other tissues (Fig. 3). 

DISCUSSION 

Many of the heat shock genes identified in this study encode 
factors that function either as molecular "chaperones" 
(HSP90a. HSP9O0. DnaJ. TCP-1) or as mediators of protein 
degradation (polvubiquitin). The identification of these se- 
quences is consistent with the biochemical basis of heat shock 
induction (10-15). Proteins undergo denaturation at elevated 
temperatures, and those that fail to maintain proper confor- 
mation must be selectively degraded (10-15). It will be inter- 
esting to determine whether the three novel heat shock- 
inducible sequences (B7-B9) mediate protein folding and 
turnover or possess some other biochemical activity. Complete 
nucleotide sequence determination, conceptual translation, 
expression monitoring, and biochemical analysis should pro^ 
vide a detailed functional understanding of these genes. 
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Fig. 3. Transcript profiles of heat shock and phorbol ester* 
regulated genes. Gene expression levels per 100.000 mRNAs U-axes) 
are shown for 15 genes (Table 1) in human bone marrow (red), brain 
(green), prostate (blue), and heart (yellow). Genes are grouped 
according to expression levels (A-C) 



Phorbol ester, a potent activator of protein kinase C ( 16. 17). 
induced a set of genes distinct from those involved in the heat 
shock pathway. The most highly induced gene identified in this 
study. /V4C-/. encodes a nuclear tyrosine kinase that may play 
a role in regulating transcription and cell cycle progression 
(18). NF-kB1, a second phorbol ester-inducible gene, is an 
intensively studied member of the Rel transcription factor 
family (19-21). The Rel proteins are activated by a large 
number of stimuli, including phorbol esters, cytokines, bacte- 
rial and viral pathogens, and ultraviolet light (19-21). Modest 
activation was observed for three sequences not known to be 
inducible by phorbol esters, including phosphoglycerate ki- 
nase, ^-microglobulin, and a novel human gene (B19). Ex- 
tensive expression monitoring with microarrays should assist in 
understanding how each of these genes integrate into the 
highly complex phorbol ester signaling pathway. 

It is striking that four novel human genes were discovered 
with an array of 1000 randomly chosen clones, particularly 
because the heat shock and phorbol ester signaling pathways 
have been so intensively studied (10-21). The facile discovery 
of these sequences underscores the fact that microarrays can 
be used for gene discovery in the absence of any sequence 
information. By this approach, clones are chosen at random 
from any library of interest and only those clones that display 
interesting expression patterns are sequenced and character- 
ized. This parallel assay, coupled with a modest DN A sequenc- 
ing facility, allows high-throughput human genome expression 
analysis and gene discovery. 

Genes that are activated or repressed by a given stimulus 
provide functional clues to the cellular pathway involved 
(22-24). Detailed examination of these gene expression "sig- 
natures" can provide a dynamic view of the mode of action of 
a given signaling substance (22-24). Microarrays may thus 
allow rapid mechanistic examination of hormones, drugs, 
elicitors. and other small molecules; moreover, functional 
analysis of transcription factors, kinases, growth factors, cyto- 
kines, receptors, and other gene products should be possible. 
Efforts are underway to develop mRNA amplification strate- 
gies to enable probe preparation from minute tissue samples. 
This capability might allow for high-throughput patient screen- 
ing in a clinical setting. 

The current detection limit of the assay allows monitoring of 
transcripts that represent -1:500.000 (wt/wt) of the total 
mRNA. This 10-fold increase in sensitivity compared with the 
original report (4) was achieved largely by modifying the 
coupling chemistry, which reduced background fluorescence. 
The significance of this improvement is considerable in that 
approximately half the human genes identified in this study, 
including all four novel sequences, exhibited expression levels 
below the original detection limit of 1:50.000 (4). 

The ability to detect 2-fold changes in expression was 
achieved by the use of two-color fluorescence in the labeling 
and detection schemes, digitized data collection, and custom 
software. The importance of this capability is underscored by 
the fact that nearly all of the genes examined here exhibited 
<6-foid changes in expression. The four novel genes, which 
showed £ 2.2-fold activation, were probably overlooked in 
previous screens that used conventional differential expression 
techniques. It may be possible to further improve the precision 
of the microarray assay by the use of closely related fluorescent 
analogs, such as Cy3 and Cy5. in the labeling and hybridization 
reactions. 

Microarrays offer a number of advantages over other po- 
tential high-capacity approaches to expression analysis. The 
chip-based approach enables small hybridization volumes, high 
array densities, and the use of fluorescence labeling and 
detection schemes. These features provide a set of pcrf r- 
mance specifications that are unattainable with filter-based 
approaches (25. 26). The use of cDNA clones provides hy- 
bridization specificity that is not readily attained with oligo- 
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nucleotide arrays (27-30). The parallel format of the assav 
provides a simultaneous differential expression readout for 
>1000 genes. This contrasts with sequencing-based methods, 
which require serial data collection for expression anarvsis (31. 
32). A commercial source of cDNA microarrays would greatly 
speed the use of a chip-based approach to expression analysis. 

The availability of large numbers of ESTs (3) provides a rich 
resource of human cDNA clones for micToarraying. The 
>400.000 ESTs in the public data bases represent a significant 
subset of all human genes (3, 33). Microarrays of thousands of 
ESTs will provide a powerful analytical tool' for future human 
gene expression studies. The -100.000 genes in the human 
genome (1 33) emphasize the need for microaiTays of greater 
density. Attempts to improve microdeposition techniques are 
underway and should allow construction of arrays containing 
a complete set of human gene targets ( http://cmgm.st anford 
edu/ -schena/). Microarrays of -100.000 cDNA elements 
would allow expression monitoring of the entire human ge- 
nome in a single hybridization. This capacity, coupled with 
detailed biochemical analysis of the individual gene products, 
would greatly speed the functional analysis of the human 
genome. 

We thank S. Ellcdge (sclledge@bcm.tmc.edu) for the human cDNA 
library, Oiagen representatives for help with plasmid purification, and 
A J. Smith and colleagues at the Protein and Nucleic Acid (PAN) 
facility (Stanford) for oligonucleotide synthesis and DNA sequencing. 
We also thank members of the Davis. Brown, and Smith laboratories 
for critical comments and helpful discussions and Synieni employees 
for technical assistance. Support for R.W.D. was provided by the 
National Science Foundation (MCB9106011) and National Institutes 
of Health (R37HG00198) and for P.O.B. bv the National Institutes of 
Health (3R21HG00450) and Howard Huehes Medical Institute. 
P.O.B. is an assistant investigator of the Howard Hughes Medical 
Institute. 

1. Watson. J. D. (1993) Gene 135. 309-315. 

2. Collins, F. S. (1995) Proc. Natl Acad. Sci. USA 92, 10821-10823. 

3. Adams, M. D.. Kelley. J. M., Gocayne, J. D., Dubnick. M., Polv- 
meropoulos. M. H., Xiao. H., Merril, C. R.. Wu, A.. Olde. B.. 
Moreno, R. F., Kerlavagc, A R.. McCombie, W. R. & Venter! 
J. C. (1991) Science 251 1651-1656. 

4. Schena. M.. Shalon. D.. Davis. R. W. & Brown, P. O (1995) 
Science 270.467-470. 

5. Shalon, D. (1996) Ph.D. thesis (Stanford University). 



Proc. Sati Acad, ScL USA 93 0996) 106 1° 

6. Schena. M. (1996) BioEuays 18, 427-431. 

7. Shalon, D.. Smith. S. J. & Brown, P. O. (1996) Genome Res. 6w 
639-645. 

8. Schena. M. Jfc Davis. R. W. (1994) Proc. Satl Acad. Sci. L'U 91. 
S39.w8397. 

9. Schena. M. & Davis. R. W. ( 1992) Proc. Sail Acad. Sci. USA 89. 
3894-3898. 

10. Jindal. S. (1996) Trends Biotechnol. 14. 17-20. 

11. Wilkinson. K. D. (1995) Annu. Rex. Nutr. IS. 161-180 

12. Jakob. U. & Buchner. J. (1994) Trends Btochem. Set. 19. 20^-211 

13. Becker. J. & Craig. E. A. (1994) Eur. J. Biochcm 219. 1 1-23. 

14. Cyr. D. M.. Unger. T. & Douglas. M. G. ( 19*4) Trends Btochem. 
Sci. 19, 176-181. 

15 365^37^ A * Weissman - J- S. & Horwich. A L. (1994) Cell 78. 

16. Newton. A C. (1995) / Bioi. Chem. 270. 28495-28498 

17. Nishizuka. Y. (1995) FASEBJ. 9. 484-496. 

18. Rohan, P. J.. Davis. P., Moskaluk. C. A. Kearns. V. . Krutzsch. 
H., Siebenlist. V. & Kelly. K. (1993) Science 259. 17 0 >-i766. 

19. Thanos. D. & Maniatis. T. (1995) Cell 80. 520-532 

20. Baeuerle. P. A & Henkel, T. (1994) Annu. Rex. Immunol. 11 
141-179. 

21. Uou. H.-C. & Baltimore. D. (1993) Curr. Optn. Celt Bioi 5. 
477-487. 

22. Cohen. G. B.. Ren. R. <fc Baltimore. D. (1995) Cell 80. 237-248. 

23. Chan. AC. Desai. D. M. & Weiss. A (1994) Annu. Rex. 
Immunol 12, 555-592. 

24. Crabtree. G. R. & Oipstone. N. A. (1994) Annu. Rex. Btochem 
63, 1045-1083. 

25. Gress, T. M., Hoheisel. J. D.. Lennon, G. G.. Zehetner, G. & 
Lehrach. H. (1992) Mamm. Genome 3, 609-619. 

26. Bernard. K.. Auphan, N.. Granjeaud, S.. Victorero. G.. Schmiti- 
Verhulst. A-M.. Jordan. B. R. & Nguven. C (1996) Nucleic Acids 
Res. 24, 1435-1442. 

27. Fodor, S. P. A. Read, J. L.. Pirrung. M. C. Strver. L. Lu. A T. 
& Solas. D. (1991) Science 251, 767-773. 

28. Southern. E. M., Maskos, U. i Elder. J. K. (1992) Genomics 13, 
1008-1017. 

29. Guo. Z., Guilfoyle, R. A. Thiel. A. J. Wang. R. & Smith. L. M. 

(1994) Nucleic Acids Res. 21 5456-5465. 

30. Matson. R. S., Rampal. J., Pentoney. S. U. Jr.. Anderson. P. D. 
A Coassin, P. (1995) Anal Biochem. 224. 110-116. 

31. Velculescu. V. E.. Zhang. L.. Vogelstein. B. &. Kimler. K.W. 

(1995) Science 270. 484-487. 

32. Adams, M. D. (1996) BtoEssavs 18. 261-262. 

33. Fields. C. Adams, M. D.. White. O. & Venter, J. C (1994) Nat. 
Genet. 7, 345-346. 







*.«• ' v- . r • 





e Genome Project adds a new dimension to Questions 
. gene expression in humans and model systems. A 
Art on page 415 summaries progr ss - the 
tenorh*bdrtis etegans Genome Project and indites 
me ways information about sequences can be used. 



Cover 




EPORTS 

tsmogenic Ages for Earthquake 447 
currence Intervals and Debris Flow Fan Depo- 
ion, Owens Valley* California 
R. Bierraan, A. R. Gillespie, M. W. Caffee 



hoautotrophic Microbial 
©systems in Deep Basalt Aquifers 
O. Stevens and J. P. McKinley 



D 4S0 



rge Arctic Temperature Change at B 455 
• Wisconsin-Holocenc Glacial Transition 

M. Cuffey, G. D. Clow, R. B. Alley, M. Scuiver, 
. D. Waddtngton, R. W. Saltus 

perplasticity in Earth's Lower Mantle: 458 
tdenee from Seismic Anisotropy and 
ck Physics 

.-L Karat , S. Zhang, H.-R. Wenk 

rge*Scale Interplanetary Magnetic 461 
Id Configurari n Revealed by Solar 
dio Bursa 

1. J. Reiner, J. Fainberg, R. G. Stone 



iASB ard f Directors I 



Newsstories, Articles, rVspectives.r^rwns.and 
Reports focus on technological developments, cfinical 
applications, and ethical concerns resuhng from the 
burgeoning of genomic information, [C efegtns tri- 
age: F. Maduro and D. Pilgrim, University of Alberta] 



Role of Yeast Insulin-Degrading Eniyme 
Homologs in Propberomone Processing and 
Bud Site Selection 

N. Adames, K. Blundell, M. N. Ashby, 

Boone 




Quantitative Monitoring of Gene 
Expression Patterns with a 
Complementary DNA Microarray 

M. Schena, D. Shalon, R. W. Davis, 

Brown 



467 



P. O. 



Gene Therapy in Peripheral Blood MK^ 470 
Lymphocytes and Bone Marrow for Vl 
ADA" lmmunodeficient Patients 

C. Bordignon, L D. Notarangelo. N. Nobiii, G. 
Ferrari, G. Casorati, P. Panina, L Maaolari, D. 
Maggioni, C Rossi, P. Servida, A. G. Ugaxio, F. 
Mavilto 

T Lymphocyte-Directed Gene 475 
Therapy for ADA" SCID: Initial U 
Trial Results After 4 Years * 1 

R. M. Blaese, K. W. Culver, A. D. Miller, C S. 
Carter, T. Fleisher, M. Clerici, G. Shearer, L 
Chang, Y. Chiang, P. Tolstoshev, J. J. Greenblart, 
S. A. .Rosenberg, H. Klein. M. Berger, C. A. 
Mullen, W. J. Ramsey, L Muul. R. A. Morgan, 
W. F. Anderson 

Physical Map and Organisation of «K 480 
Arobidopsis thaliann Chromosome 4 U 
R. Schmidt, J. West, K. Love. 11 1 

Z. Lenehan, C. Lister, H. Thompson. 

D. Boucher, C Dean 

Serial Analysis of Gene Expression /fi\ 484 
V. E. Velculescu. L Zhang, U 
B. Vogelstein. K. W. Kiruler " 

TECHNICAL COMMENTS BBSBaB 

The Radius of Gyration of an 487 
Apomyogiobin Folding Intermediate 

D. Ebeser, P. A. Jennings, P. E Wright, S. Doniach, 

K. O. Hodgson, H. Tsuruta 



D Indicates accompanying feature 



397 

Good things in 
small genomes 



7 " 



tttnng Pmtxtom 



Anna C. Rooaavaft 



*A. Latin 



JaanE. Taylor 
CrtaneHin T«n 
Nancy S. Wariar 



T tm tu rw 
Rchvri&Ncholaon 



■ sciince (issN ooss-tors) t» . MRry ori a»c*pt 

V* toat «NI In OMMttef, fry t*# Amtitcan AueciaUDtt tor ma Ad- 
wncamanl ef SctMN, 1333 H Stmai KW. MtMnpton. DC 20005. Sac 
«"<><^«Po«mo» (pufcfccoonNo. ttwtt] paid at Wasranpon. DC. and addV 
• on- *"»*ng *feoaa. C op / i yn C 1995 by 9m Anmn fturciauui tor *a 
^^ ^ f w tf Soanc a. Tha tat SCIENCE a a ragonarad fraoamam of 9m 
AMS. Do ma apc Mtma .nl raarntmha) and afnaui (51 mtuni: S97 (150 
*Banadto»*aav«on). Daman rmjamsuoachpooft (51 iwml: $22$. 
foraajn poctaoa a«r»: Manes. Carlbbatn (aurtaoa mail) $53; othar eoun- 
(air ami oaWary) $93. Ftni eau. airmail. KuMm and amarmr* ratal 
— ' " n-rtrnTli^aaar^nwiiaiiiMaLGJTiilliiJij^ 




uaa pwioi of 9m Copynpni Act to grarwad by AAAS to aaranaa 
* **** _ ^^ V *?* 0 T ?* *** c °Pmpftl Ca t a ra n c a Cam- (CCC) 7n 
ftaporvng Sarmoa. provtdad thai $3.00 pat and* la parid d-aea>r 
C<rya>jStnmS<iaraU»0197D.Th»k^< ^M ,. ^^ t^7f^ 
S075yi3 13.00. S c ene ■ t» mdasad in tna Aaooorj OuWa to 



»CCC.27 

Uh 



SCIENCE • VOL. 270 • 20 OCTOBER i995 



355 



Adlp seouence fdowng Ser** and occurs witnin 
. the ooman of Axilp thai sro*«homotogy ^ NDE 
(74). To oetete the complete S7E23 secwence and 
create the sia23A:;URA3 mutation, polymerase cnan 
reactcn (PCR) pnmen (5 '-TCGGAAGACCTCAT. 
TCTTGCTCATT7TGATATTGCTC- TGTAGATTG* 
TACTGAGAGTGCAC-3': and 5'-GCTACAAACAGC- 
GTCGACTTGAATGCCCCGACATCTTCGACTGT. 
GCGGT ATTTCAGACCG-3*) were used to amplrty 
the URA3 seouence of pRS316. and the reacDon 
cyocuawas irarg ! cm gdimo>»asttycne^eo9ane 
rectecemem |R Rothsteia Wetnods Enrymol. 1W, 
281 (1 991 )]. To create the ax/7 ^i£U? mutaten corv 
taneo on piu. a 5.0-*> Sai l fragment from oAttf 
was cioned into PUC19. and an interna* 4. f>Kb Hpa 
MCio I fragment was replaced wrtn a UBJ2 fragment 
To consiruci the ste23t::l£U2 alete la dewon cor- 
r esp ondng to 931 amno aods) cameo on p!53. a 
UEU2 fragment was used to replace the 2.8-tt> Prrt 
I-ECJ136 B tragment of STE23, which occurs within a 
6.2-kb Hind Ul-BoJ D genomic tragment earned on 
PSP72 (Promega). To create YEpMFAT, a l.6-*> 
Bam Hi tragment oomanng MFAh from pKK16 |K. 
Kuctter. R E. Sterne. J. Thomer. EMBO J. 8. 3973 
(1 989)]. was fcgated nto the Bam HI site of YEo35l p. 
E. Hffl. A. M. Myers. T. J. Koemer. A. TzagoJofl. Veasr 
2.163(1986)]. 

4. J. Chant and I. Herskowitz. Caff 65. 1203 (1991). 
:5. B. W. Matthews. Acc Cham. Res. 21.333 (1988). 

5. K. Kucfuer, H. G. Dortrnax J. Thomer. J. Caff Biol. 
120, 1203 (1993): R. Koliing and C. P. Hdtenberg, 
EMBO J. 13. 3261 (1994): C. Bertcower, D. Loayza, 
S. Michaefas. Mo/. Biol. CeB 5, 1 185 (1994). 

7. A. Benoer and J. R Pnngte. Proc. Natl. Acad. ScL 
USA 86. 9976 (1989); J. Chant. K. Corrado. J. R 
Pnngte. I. Herskowitz. Caff 65. 1213 (1991); S. 
Powers. E. Gonzaies. T. Christensen. J. Cubert. D. 
Broek, tfvd.. o. 1225: H. O. Parte J. Chant. I. Her- 
skowitz. Mature 365. 269 (1993): J. Chant. Trenos 

Genet. 10. 328 (1994): and J. R Pringte. J. 

Cefl&ol. 129. 751 (1995); J. Chant. M. Mischke. E. 
Mrtchetl. I. Herskowitz. J. R. Pnngle. ibid., p. 767. 

6. G. F. Sprague Jr.. Methods. Enzymoi. 194. 77 
(1991). 

9. Singte-tetter abbreviations for the amino acid res*- 
oues are as fouows: A. Ala; C. Cys: D. Asp: E. Gkr. F, 
Phe; G. &r. H, His; I. He: K. Lys; L Leu: M. Met; N, 
Asn; P. Pro; Q. Gin; R Arg; S. Sen T. Thr. V. Vat; W, 
Trp; and Y. Tyr. 

0. A W303 1A cwatfwe. SY2625 (M47e t/a3-J teu?-3. 
772 trpM aoe2-i canh100ssn*mta2*::FUSuiacZ 
r*s3A::FUS J -HJS3) . was the parent strain tcr the mutant 
search- SY2625 derrvatrves tor the mating assays, se- 
creted pnerumone assays, and the putse-chase exper- 
iments included the fottowng strans: Y49 lsta22-7), 
Y115 (mteli.vt£U2). Y142 &JU.1JRA3I. Y173 
(ax/7 *::LBJZl. Y220 larfJ "URA3 STe23*;UrW). Y221 
\ste23&VRA3). Y231 t*xS1&:L£U2 sTa23A.*lH/2). 
and Y233 lste23*l£U2). MAT a derivatives of 
SY2625 included the following strains: Y199 
(SY2625 made MATa), Y276 tstfi22-7). Y195 
{mtalb::LBJ2l t Y196 «fltfT*:l£U2). and Y197 
<ax/7;;UfiA3). The EG 123 (MA 7a feu? ua3 ffp7 can7 
/us4) genetic background was used to create a set of 
strains for analysis of bud site selection. EG123 de- 
rivatives included the following strains: Y175 
taxHteLSJZ). Y223 W7::UA43). Y234 (ste23A.v 
LEUZ). and Y272 {axJ1t::L£U2 st*23L':L£UZl. 
MAT* derivatives of EG123 included the following 
strains: Y214 (EG123 made MATa) and Y293 
(ax/7^.vL£U2). All strains were generated by means 
ol standard genetic or molecular methods invotang 
the appropriate constructs (23). tn particular, the ax/7 
ste23 double mutant strains were created by cross- 
ing of the appropriate MATa ste23 and MA To ax/7 
mutants, toflowed by spoliation of the resurtani Op- 
ioid and isolation of the double mutant from nonoa- 
reraaf di-fype tetrads. Gene disruptions were con- 
firmed wrth either PCR or Southern (DMA) anatysa. 

1. pl29 is 8 YEP352 |J. E. hBt. A. M. Myers. T. J. Ko- 
emer. A. Tzagotofl. Yeast Z 163(1966)) ptesmidoon- 
tamg a 5^-kb Sal I tragment of p4X17. pl5i was 
denved from pl29 by insertion of a linker at the BgJ I 
site withn AX. 7 . which led to an irvframe insertion of 
thehemaggiutnin(HA)eoitc<*(DOyp 

between amno acids 854 and 855 of the AW. 7 prod- 



uacC225isaKS^ (Stratagene) ptasmd cornering, 
a 0.5-to Bam M-Sst I fragment from pAXL 7 . Suosfc- 
tuoon mutations of the proposed acowe srte of Axflp 
were created with the use of pC225 and sne-soedfc 
rruiagenesas nvoMng aoproonaie synthetic oigenu- 
deotoes &61-H68A. 5'-GTGCTCACAAAGCGCT- 
GCCAAACCGGC-3'; ax/7-£77A, 5'-AAGAATCAT- 
GTGCGCACAAAGGTGCGC-3'; anaajtfr-E77D. 5'- 
AAGAATCATGTGATCACAAAGGTGCGC-3'). The 
mutations were cm JurneU by seouenoe anaryss. Af- 
ter mutagenese. the 0.4-to Bam Ht-Msc I tragment 
from me mutagenceo pC225 piasmids was trans- 
terred nto &XL 7 to create a set of pRS316 rjtasmos 
carryng effterent AXL1 aftetes. p!24 (axf7*H6S4). 
P130 (ax/7-£77A). and pl32 {axH-E7lOi. Sffnterty. a 
set of HA-taggeo afleies earned on YE0352 were cre- 
ated after reotacement of the p15i Bam Ht-Mse I 
fragment to generate pi6l &dl-E71A). pi 62 (ax/7- 



H6aA). and P163 Uxf7-£770i. 
32. We mank J. Becker and S. Mcnaeiis tor on^on; 
a*factor antocoes: S Mcnaens tor oiscuss:nc U r . 
putoksheo results and netpmg with the puse-cnase 
cxp enme na ; J. Brown. J. Chant, and S Sanoers 
for thee nput ooncerncng puo site selection ew ■ 
*nents: M. Raymond. F. Tamrxx and M. Whnewa\ 
for piasmias: M. Mam tor providng the STE23 
genome fragment; and H. Bussey. j. Brown. 
N. Davis. T. Favero. C de Hoog. and S. Kim tcx 
com m en ts on the manuscript. Supported by a 
grant to C.B. from me Natural Sciences and Engi- 
neering Research Council of Canada. Support tor 
M.N A was from a Cawomia Tobacco-Related Dis- 
ease Research Program postooctoral telowship 
(4FT-0083). 
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Quantitative Monitoring of Gene Expression 
Patterns with aJCompIementary DNA Microarray 

Mark Schena.* Dari Shalon/t Ronald W. Davis, 
Patrick O. Brown* 

A high-capacity system was developed to monitor the expression of many genes in 
parallel. Microarrays prepared by high-speed robotic printing of complementary DMAs n 
glass were used for quantitative expression measurements of the corresponding genes. 
Because of the small format and high density of the arrays, hybridization volumes of 2 
microliters could be used that enabled detection of rare transcripts in probe mixtures 
derived from 2 micrograms of total cellular messenger RNA. Differential expression 
measurements of 45 Arabidopsis genes were made by means of simultaneous, two-color 
fluorescence hybridization. 



temporal, developmental, topographi- 
cal, histological, and physiological patterns 
in which a gene is expressed provide clues to 
its biological role. The large and expanding 
database of complementary DNA (cDNA) 
sequences from many organisms ( 1 ) presents 
the opportunity of defining these patterns at 
the level of the whole genome. 

For these studies, we used the small flow- 
ering plant Arabidopsis xhaliana as a model 
organism. Arabidopsis possesses many ad- 
vantages for gene expression analysis, in- 
cluding the fact that it has the smallest 
genome of any higher eukaryote examined 
to date (2). Forty-five cloned Arabidopsis 
cDNAs (Table 1), including 14 complete 
sequences and 31 expressed sequence tags 
(ESTs), were used as gene-specific targets. 
We obtained the ESTs by selecting cDNA 
clones at random from an Arabidopsis 
cDNA library. Sequence analysis revealed 
that 28 of the 31 ESTs matched sequences 
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in the database (Table 1 ). Three additional 
cDNAs from other organisms served as con* 
trols in the experiments. 

The 48 cDNAs, averaging -1.0 kb, 
were amplified with the polymerase chain 
reaction (PCR) and deposited into indi- 
vidual wells of a 96-well microliter plate- 
Each sample was duplicated in two adja- 
cent wells to allow the reproducibility of 
the arraying and hybridization process t 
be tested. Samples from the microtiter 
plate were printed onto glass microscope 
slides in an area measuring 3.5 mm by 5.5 
mm with the use of a high-speed arTaying 
machine (3). The arrays were processed by 
chemical and heat treatment to attach the 
DNA sequences to the glass surface and 
denature them (3). Three arrays, printed 
in a single lot, were used for the experi- 
ments here. A single microtiter plate of 
PCR products provides sufficient material 
to prim at least 500 arrays. 

Flu rescent probes were prepared from 
total Arabidopsis mRNA (4) by a single 
round of reverse transcription (5). The Aro- 
bidopsis mRNA was supplemented with hu- 
man acetylcholine receptor ( AChR) mRNA 
at a dilution of 1 : 10,000 ( w/w) before cDNA 
synthesis, to provide an internal standard for 
calibration (5). The resulting fluorescently 
labeled cDNA mixture was hybridized to an 
array at high stringency (6) and scanned 
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h a laser (3). A high-sensitivity scan gave 
nals chat saturated the detector at nearly 
of the AToWoptsis target sites (Fig. 1A). 
tlibration relative to the AChR mRNA 
ndard (Fig. 1A) established a sensitivity 
lit of - 1 : 50,000. No detectable hybridiia- 
n was observed to either the rat glucocor- 
id receptor (Fig. 1A) or the yeast TRP4 
ig. 1A) targets even at the highest scan- 
ig sensitivity. A moderate-sensitivity scan 



of the same array allowed linear detection of 
the more abundant transcripts (Fig. IB). 
Quantitation of both scans revealed a range 
of expression levels spanning three oiden of 
magnitude for the 45 genes tested (Table 2). 
RNA blots (7) for several genes (Fig. 2) 
corroborated the expression levels measured 
with the microanay to within a factor of 5 
(Table 2). 

Differential gene expression was investi- 








1 . Gene expression monitored with the use of cDNA microarrays. Fluorescent scans represented in 
jdocotor correspond to hybndaation intensities. Color bars were calibrated from the signal obtained 
the use of known concentrations of human AChR mRNA in independent experiments. Numbers and 
rs on the axes mark the position of each cDNA. (A) High-sensitivity fluorescein scan after hybridization 
fluorescein-labeied cONA derived from wild-Type plants. (B) Same array as in (A) but scanned at 
rate sensitivity. (C and D) A single array was probed with a 1 : 1 mixture of ftuorescein-labe*ed cDNA 
wild-type plants and fissarnine-labeied cDNA from HAT4-transgenic plants. The singie array was 
scanned successively to detect the fluorescein fluorescence corresponding to mRNA from wild-type 
s (C) and the hssamme fluorescence corresponding to mRNA from HAT4-transgeruc plants (D) (E 
F) A single array was probed with a 1:1 mixture of Auorescem-labeted cDNA from root tissue and 
-nine-labeled cDNA from leaf tissue. The single array was then scanned successively to detect the 
5scein fluorescence corresponding to mRNAs expressed in roots (E) and the fcsarrune fluorescence 
■spending to mRNAs expressed in leaves (F). 



gated with a simultaneous, two-color hv 
bridization scheme, which served to mm]* 
nme experimental variation inherent in the 
comparison of independent hybriduanons 
fluorescent probes were prepared from rw 0 
mRNA sources with the use of reverse tran- 
scriptase in the presence of fluorescein- and 
lissamine-labeled nucleotide analogs, re. 
spectively (5). The two probes were then 
mixed together in equal proportions, hy- 
bridized to a single array, and scanned sep- 
arately for fluorescein and lissamine emis- 
sion after independent excitation of the two 
fluorophores (3). 

To test whether overexpression of a sin- 
gle gene could be detected in a pool of total 
Arabidopsis mRNA, we used a microaiTay to 
analy2e a transgenic line overexpressmg the 
single transcription factor HAT4 (8). Fluo- 
rescent probes representing mRNA from 
wild-type and HAT4-transgenic plants were 
labeled with fluorescein and lissamine, re- 
spectively; the two probes were then mixed 
and hybridized -to a single array. An intense 
hybridization signal was observed at the 
position of the HAT4 cDNA in the lissa- 
mine-specific scan (Fig. ID), but not in the 
fluorescein-specific scan of the same array 
(Fig. 1C). Calibration with AChR mRNA 
added to the fluorescein and lissamine 
cDNA synthesis reactions at dilutions of 
1:10,000 (Fig. 1C) and 1:100 (Fig. ID), 
respectively, revealed a 50-fold elevation of 
HAT4 mRNA in the transgenic line rela- 
tive to its abundance in wild-type piano 
(Table 2). This magnitude of HAT4 over- 
expression matched that inferred from the 
Northern (RNA) analysis within a factor of 
2 (Fig. 2 and Table 2). Expression of all the 
other genes monitored on the array differed 
by less than a factor of 5 between HAT4- 
transgenic and wild-type plants (Fig 1, C 



Wild type HATtftonsQtnfc 




1.0 0.1 0.01 ijo ai aoi 



mRNA y*g) 



Human 
AChR 




20 2J0 0.2 
mRNA (ng) 



Fig. 2. Gene expression monitored with RNA 
(Northern) blot analysis. Designated amounts of 
mRNA from wild- type and H4T4-transgeric 
plants were spotted onto nylon membranes and 
probed with the cDNAs indicated. Purified human 
AChR mRNA was used for caHjration. 



SCIENCE • VOL 270 • 20 OCTOBER 1995 



1 



and D. and Tabic 2). Hybridization of flu- 
orescein-labeled glucocorticoid receptor 
cDNA (Fig. 1C) and lissamine-labeled 
TRP4 cDNA (Fig. ID) verified the pres- 
ence of the negative control targets and the 
lack of optical cross talk between the two 
fluorophores. 

To explore a more complex alteration in 
expression patterns, we performed a second 
two-color hybridization experiment with 
fluorescein- and lissamine-labeled probes 
prepared from root and leaf mRNA, respec- 
tively. The scanning sensitivities for the 
two fluorophores were normalized by 
matching the signals resulting from AChR 



mRNA, which was added to both cDNA 
synthesis reactions at a dilution of 1:1000 
(Fig. 1, E and F). A comparison of the scans 
revealed widespread differences in gene ex- 
pression between root and leaf tissue (Fig. 1, 
E and F). The mRNA from the light-regu- 
lated CAB] gene was -500-fold more abun- 
dant in leaf (Fig. IF) than in root tissue 
(Fig. IE). The expression of 26 other genes 
differed between root and leaf tissue by 
more than a factor of 5 (Fig. 1, E and F). 

The HAT4-transgenic line we examined 
has elongated hypocotyls, early flowering, 
poor germination, and altered pigmentation 
(8). Although changes in expression were 



Table 1. Sequences contained on the cDNA microarTay. Shown is the position, the known or putative 
function, and the accession number of each cOMA in the microarray {Rg. 1 ). AH but three of the ESTs used 
in this study matched a sequence in the database. NADH. reduced form of nicotinamide adenine 
dinucieotide; ATPase. 1 adenosine triphosphatase; GTP. guanosine triphosphate. 



"osuion 


cDNA 


al.2 


AChR 


a3.4 


EST3 


a5.6 


EST6 


a7.6 


AAC1 


a9, 10 


EST12 


ail. 12 


EST13 


M.2 


CAB! 


b3. 4 


EST17 


b5. 6 


GA4 


b7, 8 


EST19 


b9. 10 


GBf-7 


d11. 12 


EST23 


c1, 2 


EST29 


c3, 4 


GBF-2 


c5. 6 


EST34 


C7, 8 


EST35 


c9, 10 


EST41 


cl 1, 12 


rGR 


dl.2 


EST42 


d3. 4 


EST45 


d5. 6 


HAT1 


d7,8 


EST46 


d9, 10 


EST49 


dll. 12 


HAT2 


el.2 


HAT 4 


e3. 4 


EST50 


e5. 6 


HATS 


7.8 


EST51 


e9, 10 


HAT22 


e11. 12 . 


EST52 


fi.2 


EST59 


f3.4 


KNATl 


15.6 


EST60 


f7.e 


EST69 


f9. 10 


PPH1 


f11, 12 


EST70 


gi.2 


EST 75 


93.4 


EST 78 


95.6 


ROC1 


97.8 


EST82 


99. 10 


EST83 


911. 12 


ESTB4 


h1.2 


EST91 


h3,4 


EST96 


h5. 6 


SARI 


h7.8 


EST100 


n9. 10 


EST103 


h1l, 12 


779P4 



Function Accession 

number 



Human AChR 
Actin 

NADH dehydrogenase 

Actin 1 

Unknown 

Actin 

Chtorophyn a/b binding 
PnosphogJycerate kinase 
GibbereUic acid biosynthesis 
Unknown 

G-box binding factor 1 
Elongation factor 
Aldolase 

G-box binding factor 2 
Chbroptast protease 
Unknown 
Catalase 

Rat glucocorticoid receptor 

Unknown 

ATPase 

Homeobox -leucine zipper 1 
Light harvesting complex 
Unknown 

Homeobox -leucine zipper 2 
Homeobox -leucine zipper 4 
Phosphoribulokinase 
Homeobox -leucine zipper 5 
Unknown 

Homeobox-teucine zipper 22 
Oxygen evolving 
Unknown 

KnorrecMike homeobox 1 
Ru&tsCO small subunrt 
Translation elongation factor 
Protein phosphatase 1 
Unknown 

ChJoropiast protease 

Unknown 

CycJophitin 

GTP binding 

Unknown 

Unknown 

Unknown 

Unknown 

Synaptobrevin 

Light harvesting complex 

Light harvesting complex 

Yeast tryptophan biosynthesis 



H36236 

227010 

M20016 

U36594t 

T45783 

M85150 

T44490 

L37126 

U36595t 

X63894 

X52256 

T04477 

X63895 

RB7034 

T14152 

T22720 

M14053 

U36596t 

J04185 

U09332 

T04063 

T 76267 

U09335 

M90394 

T04344 

M 904 16 

233675 

U09336 

T21749 

234607 

U14174 

X14564 

T42799 

U34803 

T44621 

T43698 

R65481 

LI 4644 

X59152 

233795 

T45278 

T13832 

R64816 

M90418 

218205 

X03909 

X04273 



observed for HAT«, large changes in ex. 
pression were not observed for anv ot the 
ther 44 genes we examined. This u a < 
somewhat surprising, particularly because 
comparative analysis of leaf and root tissue 
identified 27 differentially expressed genes. 
Analysts of an expanded set of genes may be 
required to identify genes whose expressi n 
changes upon HAT4 overexpression; alter- 
natively, a comparison of mRNA popula- 
tions from specific tissues of wild-type and 
HAT4-transgenic plants may allow identi- 
fication of downstream genes. 

At the current density of robotic pnnting, 
it is feasible to scale up the fabrication pro- 
cess to produce arrays containing 20,000 
cDNA targets. At this density, a single array 
would be sufficient to provide gene-specific 
targets encompassing nearly the entire rep- 
ertoire of expressed genes in the Arabidopsis 
genome (2). The availability of 20,274 ESTs 
from Arabidopsa (1,9) would provide a rich 
source of templates for such studies. 

The estimated 100,000 genes in the hu- 
man genome (JO) exceeds the number of 
Arabidopsis genes by a factor of 5 (2). This 
modest increase in complexity suggests that 
similar cDNA microarrays, prepared from 
the rapidly growing repertoire of human 
ESTs (I), could be used to determine the 
expression patterns of tens of thousands of 
human genes in diveree cell types. Coupling 
an amplification strategy to the reverse 
transcription reaction (J I) could make it 
feasible to monitor expression even in 
minute tissue samples. A wide variety of 
acute and chronic physiological and patho- 
logical conditions might lead to character* 
istic changes in the patterns of gene expres- 
sion in peripheral blood cells or other easily 
sampled tissues. In concert with cDNA mi* 
croarrays for monitoring complex expres- 
sion patterns, these tissues might therefore 
serve as sensitive in vivo sensors for clinical 
diagnosis. Microarrays of cDNAs could thus 
provide a useful link between human gene 
sequences and clinical medicine. 



Table 2. Gene expression rronrtonng by microar- 
ray and RNA Wot analyses; tg. H4T4-transgenjc 
See Table 1 for additional gene information. Ex- 
pression levels (w/w) were calibrated with the use 
tfkTK>wnarTK>umsofhurT^ 
for the microarray were determined from rrtcroar- 
ray scans Ftg. 1); values for the RNA blot were 
determined from RNA blots (Rg. 2). 



•Proprietary sequence ot Stratagene (La Jofia. CaWoma). T No mate* n the database; novel EST. 
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RNA blot 
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1:48 
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CAB! (tg) 


1:120 


1:150 


HAT4 


1:8300 


1:6300 


HAT4 (tg) 


1:150 


1210 


ROC1 


1:1200 


1:1800 


ROC1 (tg) 


1:260 


1:1300 
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Gene Therapy in Peripheral Blood 
Lymphocytes and Bone Marrow for 
ADA Immunodeficient Patients 

Claudio Bordignon/ Luigi D. Notarangelo, Nadia Nobiii, 
Giuhana. Ferrari, Giulia Casorati, Paola Panina, Evelina Mazzolari, 
Daniela Maggioni, Claudia Rossi, Paolo Servida, 
Alberto G. Ugazio, Fulvio Mavilio 

Adenosine deaminase (ADA) deficiency results in severe combined immunodeficiency 
the first genetic disorder treated by gene therapy. Two different retroviral vectors were 
used to transfer ex vivo the human ADA minigene into bone marrow cells and peripheral 
blood lymphocytes from two patients undergoing exogenous enzyme replacement ther- 
apy. After 2 years of treatment, long-term survival of T and B lymphocytes, marrow cells 
and granulocytes expressing the transferred ADA gene was demonstrated and resulted 
in normalization of the immune repertoire and restoration of cellular and humoral immunity 
After discontinuation of treatment T lymphocytes, derived from transduced peripheral 
blood lymphocytes.- were progressively replaced by marrow-derived T cells in both pa- 
tients. These results indicate successful gene transfer into long-lasting progenitor cells 
producing a functional multilineage progeny. 



Severe combined immunodeficiency asso- 
ciated with inherited deficiency of ADA 

(1) is usually fatal unless affected children 
are kept in protective isolation or the im- 
mune system is reconstituted by bone mar- 
row transplantation from a human leuko- 
cyte antigen (HLA)-identical sibling donor 

(2) . This is the therapy of choice, although 
it is available only for a minority of patients. 
In recent years, other forms of therapy have 
been developed, including transplants from 
haploidentical donors (3,4), exogenous en- 
zyme replacement (5), and somatic-cell 
gene therapy (6-9). 

We previously reported a preclinical mod- 
el in which ADA gene transfer and expression 

C. Bordignon. N. Not*, G. Ferrari. D. Maggon. C Rossi. 
P. Servida. F. Maviho. Tetetnon Gene Therapy Program 
for Genetic Diseases. DiBTT. tetnuto Soerrufco H. S. Raf- 
taete. Mfian, rtafy. 

L 0. Notarangelo. E. Mazzolari. A. G. Ugazio. Oepan- 
mem of Pediatrics, University of Bresoa MedcaJ Scnoot 
Brescia. ftaJy. 

G. Casorati. Unrta di trmmciirnce, OBIT. istrtutoSci- 
entifico H. S. RafTaeie. Mian, rtary. 
P. Parina. Roche Mtono Ricerche. Mian. ttaJy. 
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successfully restored immune functions in hu- 
man ADA-deficient (ADA") peripheral 
blood lymphocytes (PBLs) in irnmunodefi- 
ciem mice in vivo (10, J J). On the basis of 
these preclinical results, the clinical applica- 
tion of gene therapy for the treatment of 
ADA" SCID (severe combined immunodefi- 
ciency disease) patients who previously failed 
exogenous enzyme replacement therapy was 
approved by our Institutional Ethical Com- 
mittees and by the Italian National Commit- 
tee for Bioethics (12). In addition to evaluat- 
ing the safety and efficacy of the gene therapy 
procedure, the aim of the study was to define 
the relative role of PBLs and hematopoietic 
stem cells in the long-term reconstitution of 
immune functions after retroviral vector-me- 
diated ADA gene transfer. F r this purpose, 
two structurally identical vect rs expressing 
the human ADA complementary DNA 
(cDNA), distinguishable by the presence of 
alternative restriction sites in a nonfunaional 
region of the viral long-terminal repeat 
(LTR), were used to transduce PBLs and bone 
marrow (BM) cells independently. This pro- 
cedure allowed identification of the origin of 
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yjftld of the invention 

5 This invention relates to a method and apparatus 

for fabricating microarrays of biological samples for 
large scale screening assays, such as arrays of DNA 
samples to be used in DNA hybridization assays for 
genetic research and diagnostic applications. 

10 
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Background of the Invention 

A variety of methods are currently available for 
making arrays of biological macromolecules , such as 

10 arrays of nucleic acid molecules or proteins. One 
method for making ordered arrays of DNA on a porous 
membrane is a "dot blot" approach. In this method, a 
vacuum manifold transfers a plurality, e.g., 96, 
aqueous samples of DNA from 3 millimeter diameter wells 

15 to a porous membrane. A common variant of this 

procedure is a "slot-blot" method in which the wells 
have highly-elongated oval shapes. 

The DNA is immobilized on the porous membrane by 
baking the membrane or exposing it to UV radiation. 

20 This is a manual procedure practical for making one 

array at a time and usually limited to 96 samples per 
array. "Dot-blot" procedures are therefore inadequate 
for applications in which many thousand samples must be 
determined. 

25 A more efficient technique employed for making 

ordered arrays of genomic fragments uses an array of 
pins dipped into the wells, e.g., the 96 wells of a 
microtitre plate, for transferring an array of samples 
to a substrate, such as a porous membrane. One array 

30 includes pins that are designed to spot a membrane in a 
staggered fashion, for creating an array of 9216 spots 
in a 22 x 22 cm area (Lehrach, et al., 1990). A 
limitation with this approach is that the volume of DNA 
spotted in each pixel of each array is highly variable. 
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In addition, the number of arrays that can be made with 
each dipping is usually quite small. 

An alternate method of creating ordered arrays of 
nucleic acid sequences is described by Pirrung, et al. 
5 (1992), and also by Fodor, et al. (1991). The method 
involves synthesizing different nucleic acid sequences 
s at different discrete regions of a support. This 
method employs elaborate synthetic schemes, and is 
generally limited to relatively short nucleic acid 

10 sample, e.g., less than 20 bases. A related method has 
been described by Southern, et al. (1992). 

Khrapko, et al. (1991) describes a method of 
making an oligonucleotide matrix by spotting DNA onto a 
thin layer of polyacrylamide. The spotting is done 

15 manually with a micropipette. 

None of the methods or devices described in the 
prior art are designed for mass fabrication of 
microarrays characterized by (i) a large number of 
micro-sized assay regions separated by a distance of 

20 50-200 microns or less, and (ii) a well-defined amount, 
typically in the picomole range, of analyte associated 
with each region of the array. 

Furthermore, current technology is directed at 
performing such assays one at a time to a single array 

25 of DNA molecules. For example, the most common method 
for performing DNA hybridizations to arrays spotted 
onto porous membrane involves sealing the membrane in a 
plastic bag (Maniatas, et al., 1989) or a rotating 
glass cylinder (Robbins Scientific) with the labeled 

30 hybridization probe inside the sealed chamber. For 
arrays made on non-porous surfaces, such as a 
microscope slide, each array is incubated with the 
labeled hybridization probe sealed under a coverslip. 
These techniques require a separate sealed chamber for 
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each array which makes the screening and handling of 
many such arrays inconvenient and time intensive. 

Abouzied, et al. (1994) describes a method of 
printing horizontal lines of antibodies on a 
5 nitrocellulose membrane and separating regions of the 
membrane with vertical stripes of a hydrophobic 
material. Each vertical stripe is then reacted with a 
different antigen and the reaction between the 
immobilized antibody and an antigen is detected using a 

10 standard ELISA colorimetric technique. Abouzied' s 
technique makes it possible to screen many one- 
dimensional arrays simultaneously on a single sheet of 
nitrocellulose. Abouzied makes the nitrocellulose 
somewhat hydrophobic using a line drawn with PAP Pen 

15 (Research Products International) . However Abouzied 
does not describe a technology that is capable of 
completely sealing the pores of the nitrocellulose. The 
pores of the nitrocellulose are still physically open 
and so the assay reagents can leak through the 

20 hydrophobic barrier during extended high temperature 
incubations or in the presence of detergents which 
makes the Abouzied technique unacceptable for DNA 
hybridization assays. 

Porous membranes with printed patterns of 

25 hydrophilic/hydrophobic regions exist for applications 
such as ordered arrays of bacteria colonies. QA Life 
Sciences (San Diego CA) makes such a membrane with a 
grid pattern printed on it. However, this membrane has 
the same disadvantage as the Abouzied technique since 

30 reagents can still flow between the gridded arrays 
making them unusable for separate DNA hybridization 
assays . 

Pall Corporation make a 96-well plate with a 
porous filter heat sealed to the bottom of the plate. 
35 These plates are capable of containing different 
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reagents in each well without cross-contamination. 
However , each well is intended to hold only one target 
element whereas the invention described here makes a 
microarray of many biomolecules in each subdivided 
5 region of the solid support. Furthermore, the 96 well 
plates are at least 1 cm thick and prevent the use of 
the device for many color imetric, fluorescent and 
radioactive detection formats which require that the 
membrane lie flat against the detection surface. The 

10 invention described here requires no further processing 
after the assay step since the barriers elements are 
shallow and do not interfere with the detection step 
thereby greatly increasing convenience. 

Hyseq Corporation has described a method of making 

15 an " array of arrays" on a non-porous solid support for 
use with their sequencing by hybridization technique. 
The method described by Hyseq involves modifying the 
chemistry of the solid support material to form a 
hydrophobic grid pattern where each subdivided region 

20 contains a microarray of biomolecules. Hyseq 's flat 
hydrophobic pattern does not make use of physical 
blocking as an additional means of preventing cross 
contamination . 

25 Bwtim»i» Y of the Invention 

The invention includes, in one aspect, a method of 
forming a microarray of analyte-assay regions on a 
solid support, where each region in the array has a 
known amount of a selected, analyte-specif ic reagent. 

30 The method involves first loading a solution of a 
selected analyte-specif ic reagent in a reagent- 
dispensing device having an elongate capillary channel 
(i) formed by spaced-apart, coextensive elongate 
members, (ii) adapted to hold a quantity of the reagent 

35 solution and (iii) having a tip region at which aqueous 
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solution in the channel forms a meniscus. The channel 
is preferably formed by a pair of spaced-apart tapered 
elements. 

The tip of the dispensing device is tapped against 
5 a solid support at a defined position on the support 

surface with an impulse effective to break the meniscus 
in the capillary channel deposit a selected volume of 
solution on the surface, preferably a selected volume 
in the range 0.01 to 100 nl. The two steps are 

10 repeated until the desired array is formed. 

The method may be practiced in forming a plurality 
of such arrays, where the solution-depositing step is 
are applied to a selected position on each of a 
plurality of solid supports at each repeat cycle. 

15 The dispensing device may be loaded with a new 

solution, by the steps of (i) dipping the capillary 
channel of the device in a wash solution, (ii) removing 
wash solution drawn into the capillary channel, and 
(iii) dipping the capillary channel into the new 

20 reagent solution. 

Also included in the invention is an automated 
apparatus for forming a microarray of analyte-assay 
regions on a plurality of solid supports, where each 
region in the array has a known amount of a selected, 

25 analyte-specific reagent. The apparatus has a holder 
for holding, at known positions, a plurality of planar 
supports, and a reagent dispensing device of the type 
described above. 

The apparatus further includes positioning 

30 structure for positioning the dispensing device at a 
selected array position with respect to a support in 
said holder, and dispensing structure for moving the 
dispensing device into tapping engagement against a 
support with a selected impulse effective to deposit a 
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selected volume on the support, .gr., a selected volume 
in the volume range 0.01 to 100 nl. 

The positioning and dispensing structures are 
controlled by a control unit in the apparatus. The 
5 unit operates to (i) place the dispensing device at a 
loading station, (ii) move the capillary channel in the 
device into a selected reagent at the loading station, 
to load the dispensing device with the reagent, and 
(iii) dispense the reagent at a defined array position 
10 on each of the supports on said holder. The unit may 
further operate, at the end of a dispensing cycle, to 
wash the dispensing device by (i) placing the 
dispensing device at a washing station, (ii) moving the 
capillary channel in the device into a wash fluid, to 
15 load the dispensing device with the fluid, and (iii) 
remove the wash fluid prior to loading the dispensing 
device with a fresh selected reagent. 

The dispensing device in the apparatus may be one 
of a plurality of such devices which are carried on the 
20 arm for dispensing different analyte assay reagents at 
selected spaced array positions. 

In another aspect, the invention includes a 
substrate with a surface having a microarray of at 
least 10 3 distinct polynucleotide or polypeptide 
25 biopolymers in a surface area of less than about 1 cm 2 . 
Each distinct biopolymer (i) is disposed at a separate, 
defined position in said array, (ii) has a length of at 
least 50 subunits, and (iii) is present in a defined 
amount between about 0.1 femtomoles and 100 nanomoles. 
30 In one embodiment, the surface is glass slide 

surface coated with a polycationic polymer, such as 
poly lysine, and the biopolymers are polynucleotides. 
In another embodiment, the substrate has a water- 
impermeable backing, a water-permeable film formed on 
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the backing, and a grid formed on the film. The grid 
is composed of intersecting water- impervious grid 
elements extending from said backing to positions 
raised above the surface of said film, and partitions 
5 the film into a plurality of water- impervious cells. A 
biopolymer array is formed within each well. 

More generally, there is provided a substrate for 
use in detecting binding of labeled polynucleotides to 
one or more of a plurality different-sequence, 

10 immobilized polynucleotides. The substrate includes, 
in one aspect, a glass support, a coating of a 
polycationic polymer, such as poly lysine, on said 
surface of the support, and an array of distinct 
polynucleotides electrostatically bound non-covalently 

15 to said coating, where each distinct biopolymer is 

disposed at a separate, defined position in a surface 
array of polynucleotides. 

In another aspect, the substrate includes a water- 
impermeable backing, a water-permeable film formed on 

20 the backing, and a grid formed on the film, where the 
grid is composed of intersecting water-impervious grid 
elements extending from the backing to positions raised 
above the surface of the film, forming a plurality of 
cells. A biopolymer array is formed within each cell. 

25 Also forming part of the invention is a method of 

detecting differential expression of each of a 
plurality of genes in a first cell type, with respect 
to expression of the same genes in a second cell type. 
In practicing the method, there is first produced 

30 fluorescent-labeled cDNA's from mRNA's isolated from 
the two cells types, where the cDNA'S from the first 
and second cells are labeled with first and second 
different fluorescent reporters. 

A mixture of the labeled cDNA's from the two cell 

35 types is added to an array of polynucleotides 



WO 95/35505 



PCT/US95/07659 



9 

representing a plurality of known genes derived from 
the two cell types, under conditions that result in 
hybridization of the cDNA's to complementary-sequence 
polynucleotides in the array. The array is then 
5 examined by fluorescence under fluorescence excitation 
conditions in which (i) polynucleotides in the array 
that are hybridized predominantly to cDNA's derived 
from one of the first and second cell types give a 
distinct first or second fluorescence emission color, 

10 respectively, and (ii) polynucleotides in the array 

that are hybridized to substantially equal numbers of 
cDNA's derived from the first and second cell types 
give a distinct combined fluorescence emission color, 
respectively. The relative expression of known genes 

15 in the two cell types can then be determined by the 
observed fluorescence emission color of each spot. 

These and other objects and features of the 
invention will become more fully apparent when the 
following detailed description of the invention is read 

20 in conjunction with the accompanying figures. 

Brief Description of the Drawings 
Fig. 1 is a side view of a reagent-dispensing 
device having a open-capillary dispensing head 
25 constructed for use in one embodiment of the invention; 

Figs. 2A-2C illustrate steps in the delivery of a 
fixed-volume bead on a hydrophobic surface employing 
the dispensing head from Fig. 1, in accordance with one 
embodiment of the method of the invention; 
30 Fig. 3 shows a portion of a two-dimensional array 

of analyte-assay regions constructed according to the 
method of the invention; 

Fig. 4 is a planar view showing components of an 
automated apparatus for forming arrays in accordance 
35 with the inventi n. 
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Fig. 5 shows a fluorescent image of an actual 20 x 
20 array of 400 f luorescently-labeled DNA samples 
immobilized on a poly-l-lysine coated slide, where the 
total area covered by the 400 element array is 16 
5 square millimeters; 

Fig. 6 is a fluorescent image of a 1.8 cm x 1.8 cm 
microarray containing lambda clones with yeast inserts, 
the fluorescent signal arising from the hybridization 
to the array with approximately half the yeast genome 

10 labeled with a green f luorophore and the other half 
with a red f luorophore; 

Fig. 7 shows the translation of the hybridization 
image of Fig. 6 into a karyotype of the yeast genome, 
where the elements of Fig. -6 microarray contain yeast 

15 DNA sequences that have been previously physically 
mapped in the yeast genome; 

Fig. 8 show a fluorescent image of a 0.5 cmx 0.5 
cm microarray of 24 cDNA clones, where the microarray 
was hybridized simultaneously with total cDNA from wild 

20 type Arabldopsis plant labeled with a green f luorophore 
and total cDNA from a transgenic Arabidopsxs plant 
labeled with a red f luorophore, and the arrow points to 
the cDNA clone representing the gene introduced into 
the transgenic AraJbidopsis plant; 

25 Fig. 9 shows a plan view of substrate having an 

array of cells formed by barrier elements in the form 
of a grid; 

Fig. 10 shows an enlarged plan view of one of the 
cells in the substrate in Fig. 9, showing an array of 
30 polynucleotide regions in the cell; 

Fig. 11 is an enlarged sectional view of the 
substrate in Fig. 9, taken along a section line in that 
figure; and 

Fig. 12 is a scann d image of a 3 cm x 3 cm 
35 nitrocellulose solid support containing four identical 
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arrays of M13 clones in each of four quadrants, where 
each quadrant was hybridized simultaneously to a 
different oligonucleotide using an open face 
hybridization method. 

5 

Detailed Description of the Invention 

I. peEinmons 

Unless indicated otherwise, the terms defined 
below have the following meanings: 

10 "Ligand" refers to one member of a ligand/ anti- 

ligand binding pair. The ligand may be, for example, 
one of the nucleic acid strands in a complementary, 
hybridized nucleic acid duplex binding pair; an 
effector molecule in an effector /receptor binding pair; 

15 or an antigen in an antigen/ antibody or 
antigen/ antibody fragment binding pair. 

"Antiligand" refers to the opposite member of a 
ligand/anti-ligand binding pair. The antiligand may be 
the other of the nucleic acid strands in a 

20 complementary, hybridized nucleic acid duplex binding 
pair; the receptor molecule in an effector /receptor 
binding pair; or an antibody or antibody fragment 
molecule in antigen/ antibody or antigen/antibody 
fragment binding pair, respectively. 

25 "Analyte" or "analyte molecule" refers to a 

molecule, typically a macromolecule , such as a 
polynucleotide or polypeptide, whose presence, amount, 
and/ or identity are to be determined. The analyte is 
one member of a ligand/anti-ligand pair. 

30 "Analyte-specific assay reagent" refers to a 

molecule effective to bind specifically to an analyte 
molecule. The reagent is the opposite member of a 
ligand/anti-ligand binding pair. 

An "array f regions on a solid support" is a 
35 linear or two-dimensional array of preferably discrete 
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regions, each having a finite area, formed on the 
surface of a solid support • 

A "microarray" is an array of regions having a 
density of discrete regions of at least about 100/cm 2 , 
5 and preferably at least about 1000 /cm 2 . The regions in 
a microarray have typical dimensions, e.g., diameters, 
in the range of between about 10-250 jim, and are 
separated from other regions in the array by about the 
same distance. 

10 A support surface is "hydrophobic" if a aqueous- 

medium droplet applied to the surface does not spread 
out substantially beyond the area size of the applied 
droplet. That is, the surface acts to prevent 
spreading of the droplet applied to the surface by 

15 hydrophobic interaction with the droplet. 

A "meniscus" means a concave or convex surface 
that forms on the bottom of a liquid in a channel as a 
result of the surface tension of the liquid. 

"Distinct biopolymers", as applied to the 

20 biopolymers forming a microarray, means an array member 
which is distinct from other array members on the basis 
of a different biopolymer sequence, and/ or different 
concentrations of the same or distinct biopolymers, 
and/or different mixtures of distinct or different- 

25 concentration biopolymers . Thus an array of "distinct 
polynucleotides" means an array containing, as its 
members, (i) distinct polynucleotides, which may have a 
defined amount in each member, (ii) different, graded 
concentrations of given-sequence polynucleotides, 

30 and/or (iii) different-composition mixtures of two or 
more distinct polynucleotides. 

"Cell type" means a cell from a given source, 
e.g., a tissue, or organ, or a cell in a given state of 



WQ 95/35505 



PCT/US95/07659 



13 

differentiation, or a cell associated with a given 
pathology or genetic makeup. 

II. Method of Microarrav Formation 
5 This section describes a method of forming a 

microarray of analyte-assay regions on a solid support 
or substrate, where each region in the array has a 
known amount of a selected, analyte-specif ic reagent. 

Fig. 1 illustrates, in a partially schematic view, 

10 a reagent-dispensing device 10 useful in practicing the 
method. The device generally includes a reagent 
dispenser 12 having an elongate open capillary channel 
14 adapted to hold a quantity of the reagent solution, 
such as indicated at 16, as will be described below. 

15 The capillary channel is formed by a pair of spaced- 

apart, coextensive, elongate members 12a, 12b which are 
tapered toward one another and converge at a tip or tip 
region 18 at the lower end of the channel. More 
generally, the open channel is formed by at least two 

20 elongate, spaced-apart members adapted to hold a 

quantity of reagent solutions and having a tip region 
at which aqueous solution in the channel f orms a 
meniscus, such as the concave meniscus illustrated at 
20 in Fig. 2A. The advantages of the open channel 

-25 construction of the dispenser are discussed below. 

With continued reference to Fig. 1, the dispenser 
device also includes structure for moving the dispenser 
rapidly toward and away from a support surface, for 
effecting deposition of a known amount of solution in 

30 the dispenser on a support, as will be described below 
with reference to Figs. 2A-2C. In the embodiment 
shown, this structure includes a solenoid 22 which is 
activatable to draw a solenoid piston 24 rapidly 
downwardly, then release the piston, .g., under spring 

35 bias, to a normal, raised position, as shown. The 
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dispenser is carried on the piston by a connecting 
member 26, as shown. The just-described moving 
structure is also referred to herein as dispensing 
means for moving the dispenser into engagement with a 
5 solid support, for dispensing a known volume of fluid 
on the support. 

The dispensing device just described is carried on 
an arm 28 that may be moved either linearly or in an x- 
y plane to position the dispenser at a selected 

10 deposition position, as will be described. 

Figs. 2A-2C illustrate the method of depositing a 
known amount of reagent solution in the just-described 
dispenser on the surface of a solid support, such as 
the support indicated at 30. The support is a polymer, 

15 glass, or other solid-material support having a surface 
indicated at 31. 

In one general embodiment, the surface is a 
relatively hydrophilic, I.e., wettable surface, such as 
a surface having native, bound or covalently attached 

20 charged groups. On such surface described below is a 
glass surface having an absorbed layer of a 
polycationic polymer, such as poly-l-lysine. 

In another embodiment, the surface has or is 
formed to have a relatively hydrophobic character, 

25 i.e., one that causes aqueous medium deposited on the 
surface to bead. A variety of known hydrophobic 
polymers, such as polystyrene, polypropylene, or 
polyethylene have desired hydrophobic properties, as do 
glass and a variety of lubricant or other hydrophobic 

30 films that may be applied, to the support surface. 

Initially, the dispenser is loaded with a selected 
analyte-specif ic reagent solution, such as by dipping 
the dispenser tip, after washing, into a solution of 
the reagent, and allowing filling by capillary flow 

35 into the dispenser channel. The dispenser is now moved 
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to a selected position with respect to a support 
surface, placing the dispenser tip directly above the 
support-surface position at which the reagent is to be 
deposited. This movement takes place with the 
5 dispenser tip in its raised position, as seen in Fig. 
2A, where the tip is typically at least several 1-5 mm 
above the surface of the substrate. 

With the dispenser so positioned, solenoid 22 is 
now activated to cause the dispenser tip to move 
10 rapidly toward and away from the substrate surface, 
making momentary contact with the surface, in effect, 
tapping the tip of the dispenser against the support 
surface. The tapping movement of the tip against the 
surface acts to break the liquid meniscus in the tip 
15 channel, bringing the liquid in the tip into contact 
with the support surface. This, in turn, produces a 
flowing of the liquid into the capillary space between 
the tip and the surface, acting to draw liquid out of 
the dispenser channel, as seen in Fig. 2B. 
20 Fig. 2C shows flow of fluid from the tip onto the 

support surface, which in this case is a hydrophobic 
surface. The figure illustrates that liquid continues 
to flow from the dispenser onto the support surface 
until it forms a liquid bead 32. At a given bead size, 
25 i.e., volume, the tendency of liquid to flow onto *the 
surface will be balanced by the hydrophobic surface 
interaction of the bead with the support surface, which 
acts to limit the total bead area on the surface, and 
by the surface tension of the droplet, which tends 
30 toward a given bead curvature. At this point, a given 
bead volume will have formed, and continued contact of 
the dispenser tip with the bead, as the dispenser tip 
is being withdrawn, will have little or no effect on 
bead volume. 
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For liquid-dispensing on a more hydrophilic 
surface, the liquid will have less of a tendency to 
bead, and the dispensed volume will be more sensitive 
to the total dwell time of the dispenser tip in the 
5 immediate vicinity of the support surface, e.gr., the 
positions illustrated in Figs. 2B and 2C. 

The desired deposition volume, i.e., bead volume, 
formed by this method is preferably in the range 2 pi 
(picoliters) to 2 nl (nanoliters) , although volumes as 

10 high as 100 nl or more may be dispensed. It will be 
appreciated that the selected dispensed volume will 
depend on (i) the "footprint" of the dispenser tip, 
i.e., the size of the area spanned by the tip, (ii) the 
hydrophobicity of the support surface, and (iii) the 

15 time of contact with and rate of withdrawal of the tip 
from the support surface. In addition, bead size may 
be reduced by increasing the viscosity of the medium, 
effectively reducing the flow time of liquid from the 
dispenser onto the support surface. The drop size may 

20 be further constrained by depositing the drop in a 
hydrophilic region surrounded by a hydrophobic grid 
pattern on the support surface. 

In a typical embodiment, the dispenser tip is 
tapped rapidly against the support surface, with a 

25 total residence time in contact with the support of 
less than about 1 msec, and a rate of upward travel 
from the surface of about 10 cm/sec. 

Assuming that the bead that forms on contact with 
the surface is a hemispherical bead, with a diameter 

30 approximately equal to the width of the dispenser tip, 
as shown in Fig. 2C, the volume of the bead formed in 
relation to dispenser tip width (d) is given in Table l 
below. As seen, the volume of the bead ranges between 
2 pi to 2 nl as the width size is increased from about 

35 20 to 200 nm. 
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Table 1 



d 


Volume (nl) 


20 fim 


2 x 10 ' 3 


50 fim 


3.1 x lO' 2 


100 urn 


2.5 x 10" 1 


200 Jim 


2 



10 At a given tip size, bead volume can be reduced in 

a controlled fashion by increasing surface 
hydrophobicity, reducing time of contact of the tip 
with the surface, increasing rate of movement of the 
tip away from the surface, and/ or increasing the 

15 viscosity of the medium. Once these parameters are 

fixed, a selected deposition volume in the desired pi 
to nl range can be achieved in a repeatable fashion. 

After depositing a bead at one selected location 
on a support, the tip is typically moved to a 

20 corresponding position on a second support, a droplet 
is deposited at that position, and this process is 
repeated until a liquid droplet of the reagent has been 
deposited at a selected position on each of a plurality 
of supports. 

25 The tip is then washed to remove the reagent 

liquid, filled with another reagent liquid and this 
reagent is now deposited at each another array position 
on each of the supports. In one embodiment, the tip is 
washed and refilled by the steps of (i) dipping the 

30 capillary channel of the device in a wash solution, 
(ii) removing wash solution drawn into the capillary 
channel, and (iii) dipping the capillary channel into 
the new reagent solution. 

From the foreg ing, it will be appreciated that 

35 the tweezers-like, open-capillary dispenser tip 
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provides the advantages that (i) the open channel of 
the tip facilitates rapid, efficient washing and drying 
before reloading the tip with a new reagent, (ii) 
passive capillary action can load the sample directly 
5 from a standard microwell plate while retaining 

sufficient sample in the open capillary reservoir for 
the printing of numerous arrays, (iii) open capillaries 
are less prone to clogging than closed capillaries, and 
(iv) open capillaries do not require a perfectly faced 

10 bottom surface for fluid delivery. 

A portion of a microarray 36 formed on the surface 
38 of a solid support 40 in accordance with the method 
just described is shown in Fig. 3. The array is formed 
of a plurality of analyte-specif ic reagent regions, 

15 such as regions 42, where each region may include a 
different analyte-specif ic reagent. As indicated 
above, the diameter of each region is preferably 
between about 20-200 /ra. The spacing between each 
region and its closest (non-diagonal) neighbor, 

20 measured from center-to-center (indicated at 44) , is 

preferably in the range of about 20-400 pa. Thus, for 
example, an array having a center-to-center spacing of 
about 250 contains about 40 regions/cm or 1,600 
regions/cm 2 . After formation of the array, the support 

25 is treated to evaporate the liquid of the droplet 

forming each region, to leave a desired array of dried, 
relatively flat regions. This drying may be done by 
heating or under vacuum. 

In some cases, it is desired to first rehydrate 

30 the droplets containing the analyte reagents to allow 
for more time for adsorption to the solid support. It 
is also possible to spot out the analyte reagents in a 
humid environment so that droplets do not dry until the 
arraying operation is complete. 
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III. Automated Apparatus for Forming Arrays 

In another aspect, the invention includes an 
automated apparatus for forming an array of analyte- 
assay regions on a solid support, where each region in 
5 the array has a known amount of a selected, analyte- 
specific reagent. 

The apparatus is shown in planar, and partially 
schematic view in Fig. 4. A dispenser device 72 in the 
apparatus has the basic construction described above 

10 with respect to Fig. 1, and includes a dispenser 74 

having an open-capillary channel terminating at a tip, 
substantially as shown in Figs. 1 and 2A-2C. 

The dispenser is mounted in the device for 
movement toward and away from a dispensing position at 

15 which the tip of the dispenser taps a support surface, 
to dispense a selected volume of reagent solution, as 
described above. This movement is effected by a 
solenoid 76 as described above. Solenoid 76 is under 
the control of a control unit 77 whose operation will 

20 be described below. The solenoid is also referred to 
herein as dispensing means for moving the device into 
tapping engagement with a support, when the device is 
positioned at a defined array position with respect to 
that support. 

25 The dispenser device is carried on an arm 74 which 

is threadedly mounted on a worm screw 80 driven 
(rotated) in a desired direction by a stepper motor 82 
also under the control of unit 77. At its left end in 
the figure screw 80 is carried in a sleeve 84 for 

30 rotation about the screw axis. At its other end, the 
screw is mounted to the drive shaft of the stepper 
motor, which in turn is carried on a sleeve 86. The 
dispenser device, worm screw, the two sleeves mounting 
the w rm screw, and the stepper motor used in moving 

35 the device in the "x w (horizontal) direction in the 
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figure form what is referred to here collectively as a 
displacement assembly 86. 

The displacement assembly is constructed to 
produce precise, micro-range movement in the direction 
5 of the screw, i.e., along an x axis in the figure. In 
one mode, the assembly functions to move the dispenser 
in x-axis increments having a selected distance in the 
range 5-25 /xm. In another mode, the dispenser unit may 
be moved in precise x-axis increments of several 

10 microns or more,; for positioning the dispenser at 

associated positions on adjacent supports, as will be 
described below. 

The displacement assembly, in turn, is mounted for 
movement in the M y" (vertical) axis of the figure, for 

15 positioning the dispenser at a selected y axis 

position. The structure mounting the assembly includes 
a fixed rod 88 mounted rigidly between a pair of frame 
bars 90, 92, and a worm screw 94 mounted for rotation 
between a pair of frame bars 96, 98. The worm screw is 

20 driven (rotated) by a stepper motor 100 which operates 
under the control of unit 77. The motor is mounted on 
bar 96, as shown. 

The structure just described, including worm screw 
94 and motor 100, is constructed to produce precise, 

25 micro-range movement in the direction of the screw, 
i.e., along an y axis in the figure. As above, the 
structure functions in one mode to move the dispenser 
in y-axis increments having a selected distance in the 
range 5-250 fim, and in a second mode, to move the 

30 dispenser in precise y-axis increments of several 

microns (pm) or more, for positioning the dispenser at 
associated positions on adjacent supports. 

The displacement assembly and structure for moving 
this assembly in the y axis are referred to herein 

35 coll ctively as positioning means for positioning the 
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dispensing device at a selected array position with 
respect to a support. 

A holder 102 in the apparatus functions to hold a 
plurality of supports, such as supports 104 on which 
5 the microarrays of regent regions are to be formed by 
the apparatus. The holder provides a number of 
recessed slots, such as slot 106, which receive the 
supports, and position them at precise selected 
positions with respect to the frame bars on which the 

10 dispenser moving means is mounted. 

As noted above, the control unit in the device 
functions to actuate the two stepper motors and 
dispenser solenoid in a sequence designed for automated 
operation of the apparatus in forming a selected 

15 microarray of reagent regions on each of a plurality of 
supports . 

The control unit is constructed, according to 
conventional microprocessor control principles, to 
provide appropriate signals to each of the solenoid and 
20 each of the stepper motors, in a given timed sequence 
and for appropriate signalling time. The construction 
of the unit, and the settings that are selected by the 
user to achieve a desired array pattern, will be 
understood from the following description of a typical 
25 apparatus operation. 

Initially, one or more supports are placed in one 
or more slots in the holder. The dispenser is then 
moved to a position directly above a well (not shown) 
containing a solution of the first reagent to be 
30 dispensed on the support (s) . The dispenser solenoid is 
actuated now to lower the dispenser tip into this well, 
causing the capillary channel in the dispenser to fill. 
Motors 82, 100 are now actuated to position the 
dispenser at a s lected array position at the first of 
35 the supports. Solenoid actuation of the dispenser is 
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then effective to dispens a selected-volume droplet of 
that reagent at this location. As noted above, this 
operation is effective to dispense a selected volume 
preferably between 2 pi and 2 nl of the reagent 
5 solution. 

The dispenser is now moved to the corresponding 
position at an adjacent support and a similar volume of 
the solution is dispensed at this position. The 
process is repeated until the reagent has been 

10 dispensed at this preselected corresponding position on 
each of the supports. 

Where it is desired to dispense a single reagent 
at more than two array positions on a support, the 
dispenser may be moved to different array positions at 

15 each support, before moving the dispenser to a new 
support, or solution can be dispensed at individual 
positions on each support, at one selected position, 
then the cycle repeated for each new array position. 
To dispense the next reagent, the dispenser is 

20 positioned over a wash solution (not shown) , and the 
dispenser tip is dipped in and out of this solution 
until the reagent solution has been substantially 
washed from the tip. Solution can be removed from the 
tip, after each dipping, by vacuum, compressed air 

25 spray, sponge, or the like. 

The dispenser tip is now dipped in a second 
reagent well, and the filled tip is moved to a second 
selected array position in the first support. The 
process of dispensing reagent at each of the 

30 corresponding second-array positions is then carried as 
above. This process is repeated until an entire 
microarray of reagent solutions on each of the supports 
has been formed. 

35 IV. MicroqSTay Substrate 
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This section describes embodiments of a substrate 
having a microarray of biological polymers carried on 
the substrate surface. Subsection A describes a multi- 
cell substrate, each cell of which contains a 
5 microarray, and preferably an identical microarray, of 
distinct biopolymers, such as distinct polynucleotides, 
formed on a porous surface. Subsection B describes a 
microarray of distinct polynucleotides bound on a glass 
slide coated with a polycationic polymer. 

10 

A. Multi-CeU Substrate 

Fig. 9 illustrates, in plan view, a substrate 110 
constructed according to the invention. The substrate 
has an 8 x 12 rectangular array 112 of cells, such as 

15 cells 114, 116, formed on the substrate surface. With 
. reference to Fig. 10, each cell, such as cell 114, in 
turn supports a microarray 118 of distinct biopolymers, 
such as polypeptides or polynucleotides at known, 
addressable regions of the microarray. Two such 

20 regions forming the microarray are indicated at 120, 

and correspond to regions, such as regions 42, forming 
the microarray of distinct biopolymers shown in Fig. 3. 

The 96-cell array shown in Fig. 9 has typically 
array dimensions between about 12 and 244 mm in width 

25 and 8 and 400 mm in length, with the cells in the array 
having width and length dimension of 1/12 and 1/8 the 
array width and length dimensions, respectively, i.e., 
between about 1 and 20 in width and 1 and 50 mm in 
length. 

30 The construction of substrate is shown cross- 

sect ionally in Fig. 11, which is an enlarged sectional 
view taken along view line 124 in Fig. 9. The 
substrate includes a water-impermeable backing 126, 
such as a glass slide or rigid polymer sheet. Formed 

35 on the surface of the backing is a water-permeable film 
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128. The film is form d of a porous membrane material, 
such as nitrocellulose membrane, or a porous web 
material, such as a nylon, polypropylene, or PVDF 
porous polymer material. The thickness of the film is 
5 preferably between about 10 and 1000 /im. The film may 
be applied to the backing by spraying or coating 
uncured material on the backing, or by applying a 
preformed membrane to the backing. The backing and 
film may be obtained as a preformed unit from 

10 commercial source, e.g., a plastic-backed 

nitrocellulose film available from Schleicher and 
Schuell Corporation. 

With continued reference to Fig. 11, the film- 
covered surface in the substrate is partitioned into a 

15 desired array of cells by water- impermeable grid lines, 
such as lines 130, 132, which have infiltrated the film 
down to the level of the backing, and extend above the 
surface of the film as shown, typically a distance of 
100 to 2000 /im above the film surface. 

20 The grid lines are formed on the substrate by 

laying down an uncured or otherwise f lowable resin or 
elastomer solution in an array grid, allowing the 
material to infiltrate the porous film down to the 
backing, then curing or otherwise hardening the grid 

25 lines to form the cell-array_substrate . 

One preferred material for the grid is a f lowable 
silicone available from Loctite Corporation. The 
barrier material can be extruded through a narrow 
syringe (e.g., 22 gauge) using air pressure or 

30 mechanical pressure. The syringe is moved relative to 
the solid support to print the barrier elements as a 
grid pattern. The extruded bead of silicone wicks into 
the pores of the solid support and cures to form a 
shallow waterproof barrier separating the regions of 

35 the solid support. 
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In alternativ embodiments, the barrier element 
can be a wax-based material or a thermoset material 
such as epoxy. The barrier material can also be a UV- 
curing polymer which is exposed to UV light after being 
5 printed onto the solid support. The barrier material 
may also be applied to the solid support using printing 
techniques such as silk-screen printing. The barrier 
material may also be a heat-seal stamping of the porous 
solid support which seals its pores and forms a water- 

10 impervious barrier element. The barrier material may 
also be a shallow grid which is laminated or otherwise 
adhered to the solid support. 

In addition to plastic-backed nitrocellulose, the 
solid support can be virtually any porous membrane with 

15 or without a non-porous backing. Such membranes are 
readily available from numerous vendors and are made 
from nylon, PVDF, polysulfone and the like. In an 
alternative embodiment, the barrier element may also be 
used to adhere the porous membrane to a non-porous 

20 backing in addition to functioning as a barrier to 
prevent cross contamination of the assay reagents. 

In an alternative embodiment, the solid support 
can be of a non-porous material. The barrier can be 
printed either before or after the microarray of 

25 biomolecules is printed on the solid support. 

As can be appreciated, the cells formed by the 
grid lines and the underlying backing are water- 
impermeable, having side barriers projecting above the 
porous film in the cells. Thus, def ined-volume samples 

30 can be placed in each well without risk of cross- 
contamination with sample material in adjacent cells. 
In Fig. 11, defined volumes samples, such as sample 
134, are shown in the cells. 

As noted above, each well contains a microarray of 

35 distinct biopolymers. In one general embodiment, th 
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microarrays in the well are identical arrays of 
distinct biopolymers, e.g., different sequence 
polynucleotides. Such arrays can be formed in 
accordance with the methods described in Section II, by 
5 depositing a first selected polynucleotide at the seme 
selected microarray position in each of the cells, then 
depositing a second polynucleotide at a different 
microarray position in each well, and so on until a 
complete, identical microarray is formed in each cell. 

10 In a preferred embodiment, each microarray 

contains about 10 3 distinct polynucleotide or 
polypeptide biopolymers per surface area of less than 
about 1 cm 2 . Also in a preferred embodiment, the 
-biopolymers in each microarray region are present in a 

15 defined amount between about 0.1 femtomoles and 100 

nanomoles. The ability to form high-density arrays of 
biopolymers, where each region is formed of a well- 
defined amount of deposited material, can be achieved 
in accordance with the microarr ay-forming method 

20 described in Section II. 

Also in a preferred embodiments, the biopolymers 
are polynucleotides having lengths of at least about 50 
bp, i.e., substantially longer than oligonucleotides 
which can be formed in high-density arrays by schemes 

25 involving parallel, step-wise polymer synthesis on the 
array surface. 

In the case of a polynucleotide array, in an assay 
procedure, a small volume of the labeled DNA probe 
mixture in a standard hybridization solution is loaded 

30 onto each cell. The solution will spread to cover the 
entire microarray and stop at the barrier elements. 
The solid support is then incubated in a humid chamber 
at the appropriate temperature as required by the 
assay. 



WO 95/35505 



PCT/US95/07659 



27 



Each assay may be conducted in an "open-face" 
format where no further sealing step is required, since 
the hybridization solution will be kept properly 
hydrated by the water vapor in the humid chamber. At 

5 the conclusion of the incubation step, the entire solid 
support containing the numerous microarrays is rinsed 
quickly enough to dilute the assay reagents so that no 
significant cross contamination occurs. The entire 
solid support is then reacted with detection reagents 
10 if needed and analyzed using standard color imetric , 
radioactive or fluorescent detection means. All . 
processing and detection steps are performed 
simultaneously to all of the microarrays on the solid 
support ensuring uniform assay conditions for all of 

15 the microarrays on the solid support. 

B. Glass-Sljde Polyn ucleotide Array 
Fig. 5 shows a substrate 136 formed according to 
another aspect of the invention, and intended for use 
20 in detecting binding of labeled polynucleotides to one 
or more of a plurality distinct polynucleotides. The 
substrate includes a glass substrate 138 having formed 
on its surface, a coating of a polycationic polymer, 
preferably a cationic polypeptide, such as polylysine 
25 or polyarginine. Formed on the polycationic coating is 
a microarray 140 of distinct polynucleotides, each 
localized at known selected array regions, such as 
regions 142. 

The slide is coated by placing a uniform-thickness 
30 film of a polycationic polymer, e.g., poly-l-lysine, on 
the surface of a slide and drying the film to form a 
dried coating. The amount of polycationic polymer 
added is sufficient to form at least a monolayer of 
polymers on the glass surface. The polymer film is 
35 bound to surfac via electrostatic binding between 
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negative silyl-OH groups on the surface and charged 
amine groups in the polymers. Poly-l-lysine coated 
glass slides may be obtained commercially, e.g., from 
Sigma Chemical Co. (St. Louis, MO). 
5 To form the microarray, defined volumes of 

distinct polynucleotides are deposited on the polymer- 
coated slide, as described in Section II. According to 
an important feature of the substrate, the deposited 
polynucleotides remain bound to the coated slide 

10 surface non-covalently when an aqueous DNA sample is 
applied to the substrate under conditions which allow 
hybridization of reporter-labeled polynucleotides in 
the sample to complementary-sequence (single-stranded) 
polynucleotides in the substrate array. The method is 

15 illustrated in Examples 1 and 2. 

To illustrate this feature, a substrate of the 
type just described, but having an array of same- 
sequence polynucleotides, was mixed with fluorescent- 
labeled complementary DNA under hybridization 

20 conditions. After washing to remove non-hybridized 
material, the substrate was examined by low-power 
fluorescence microscopy. The array can be visualized 
by the relatively uniform labeling pattern of the array 
regions. 

25 in a preferred embodiment, each microarray 

contains at least 10 3 distinct polynucleotide or 
polypeptide biopolymers per surface area of less than 
about 1 cm 2 . In the embodiment shown in Fig. 5, the 
microarray contains 400 regions in an area of about 16 

30 mm 2 , or 2.5 x 10 3 regions/cm 2 . Also in a preferred 

embodiment, the polynucleotides in the each microarray 
region are present in a defined amount between about 
0.1 femtomoles and 100 nanomoles in the case of 
polynucleotides. As above, the ability to form high- 
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density arrays of this type, where each region is 
formed of a well-defined amount of deposited material, 
can be achieved in accordance with the microarray- 
forming method described in Section II. 
5 Also in a preferred embodiments, the 

polynucleotides have lengths of at least about 50 bp, 
i.e., substantially longer than oligonucleotides which 
can be formed in high-density arrays by various in situ 
synthesis schemes. 

10 

V. Utility 

Microarrays of immobilized nucleic acid sequences 
prepared in accordance with the invention can be used 
for large scale hybridization assays in numerous 

15 genetic applications, including genetic and physical 

mapping of genomes, monitoring of gene expression, DNA 
sequencing, genetic diagnosis, g en o typing of organisms, 
and distribution of DNA reagents to researchers. 

For gene mapping, a gene or a cloned DNA fragment 

20 is hybridized to an ordered array of DNA fragments, and 
the identity of the DNA elements applied to the array 
is unambiguously established by the pixel or pattern of 
pixels of the array that are detected. One application 
of such arrays for creating a genetic map is described 

25 by Nelson, et al. (1993). In constructing physical 
maps of the genome, arrays of immobilized cloned DNA 
fragments are hybridized with other cloned DNA 
fragments to establish whether the cloned fragments in 
the probe mixture overlap and are therefore contiguous 

30 to the immobilized clones on the array. For example, 
Lehrach, et al., describe such a process. 

The arrays of immobilized DNA fragments may also 
be used for genetic diagnostics. To illustrate, an 
array containing multiple forms f a mutated gene or 

35 genes can be pr bed with a labeled mixture fa 
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patient's DNA which will preferentially interact with 
only one of the immobilized versions of the gene. 

The detection of this interaction can lead to a 
medical diagnosis. Arrays of immobilized DNA fragments 
5 can also be used in DNA probe diagnostics. For 

example, the identity of a pathogenic microorganism can 
be established unambiguously by hybridizing a sample of 
the unknown pathogen's DNA to an array containing many 
types of known pathogenic DNA. A similar technique can 

10 also be used for .unambiguous genotyping of any 

organism. Other molecules of genetic interest, such as 
cDNA's and PNA's can be immobilized on the array or 
alternately used as the labeled probe mixture that is 
applied to the array. 

15 In one application, an array of cDNA clones 

representing genes is hybridized with total cDNA from 
an organism to monitor gene expression for research or 
diagnostic purposes. Labeling total cDNA from a normal 
cell with one color f luorophore and total cDNA from a 

20 diseased cell with another color f luorophore and 

simultaneously hybridizing the two cDNA samples to the 
same array of cDNA clones allows for differential gene 
expression to be measured as the ratio of the two 
f luorophore intensities. This two-color experiment can 

25 be used to monitor gene expression in different tissue 
types, disease states, response to drugs, or response 
to environmental factors. & An example of this approach 
is illustrated in Examples 2 , described with respect to 
Fig. 8. 

30 By way of example and without implying a 

limitation of scope, such a procedure could be used to 
simultaneously screen many patients against all known 
mutations in a disease gene. This invention could be 
used in the form of, for example, 96 identical 0.9 cm x 

35 2.2 cm microarrays fabricated on a single 12 cm x 18 cm 
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sheet of plastic-backed nitrocellulose where each 
microarray could contain, for example, 100 DNA 
fragments representing all known mutations of a given 
gene. The region of interest from each of the DNA 
5 samples from 96 patients could be amplified, labeled, 
and hybridized to the 96 individual arrays with each 
assay performed in 100 microliters of hybridization 
solution. The approximately 1 thick silicone rubber 
barrier elements between individual arrays prevent 

10 cross contamination of the patient samples by sealing 
the pores of the nitrocellulose and by acting as a 
physical barrier between each microarray. The solid 
support containing all 96 microarrays assayed with the 
96 patient samples is incubated, rinsed, detected and 

15 analyzed as a single sheet of material using standard 
radioactive, fluorescent, or colorimetric detection 
means (Maniatas, et al., 1989). Previously, such a 
procedure would involve the handling, processing and 
tracking of 96 separate membranes in 96 separate sealed 

20 chambers. By processing all 96 arrays as a single 

sheet of material, significant time and cost savings 
are possible. 

The assay format can be reversed where the patient 
or organism's DNA is immobilized as the array elements 

25 and each array is hybridized with a different mutated 
allele or genetic marker. The gridded solid support 
can also be used for parallel non-DNA ELISA assays. 
Furthermore, the invention allows for the use of all 
standard detection methods without the need to remove 

30 the shallow barrier elements to carry out the detection 
step. 

In addition to the genetic applications listed 
above, arrays of whole cells, peptides, enzymes, 
antibodies, antigens, receptors, ligands, 
35 phospholipids, polymers, drug cogener preparations or 
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described in this invention for large scale screening 
assays in medical diagnostics, drug discovery, 
molecular biology, immunology and toxicology, 
5 The multi-cell substrate aspect of the invention 

allows for the rapid and convenient screening of many 
DNA probes against many ordered arrays of DNA 
fragments. This eliminates the need to handle and 
detect many individual arrays for performing mass 
10 screenings for genetic research and diagnostic 

applications. Numerous microarrays can be fabricated 
on the same solid support and each microarray reacted 
with a different DNA probe while the solid support is 
processed as a single sheet of material. 

15 

The following examples illustrate, but in no way 
are intended to limit, the present invention.* 



Example 1 

20 Genomic-Complexitv Hybridization to Micro 

DNA Arrays Representing the Yeast 
Saccharomyces cerevisiae Genome with 
Two-Color Fluorescent Detection 

The array elements were randomly amplified PCR 

25 (Bohlander, et al., 1992) products using physically 

mapped lambda clones of S, cerevisiae genomic DNA 

templates (Riles, et al., 1993). The PCR was performed 

directly on the lambda phage lysates resulting in an 

amplification of both the 35 kb lambda vector and the 

30 5-15 kb yeast insert sequences in the form of a uniform 

distribution of PCR product between 250-1500 base pairs 

in length. The PCR product was purified using 

Sephadex G50 gel filtration (Pharmacia, Piscataway, NJ) 

and concentrated by evaporation to dryness at room 

35 temperature overnight. Each of the 864 amplified 
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lambda clones was rehydrated in 15 pi of 3 * SSC in 
preparation for spotting onto the glass. 

The micro arrays were fabricated on microscope 
slides which were coated with a layer of poly-l-lysine 
5 (Sigma) . The automated apparatus described in Section 
IV loaded 1 pi of the concentrated lambda clone PCR 
product in 3 x SSC directly from 96 well storage plates 
into the open capillary printing element and deposited 
-5 nl of sample per slide at 380 micron spacing between 
10 spots, on each of 40 slides* The process was repeated 
for all 864 samples and 8 control spots. After the 
spotting operation was complete, the slides were 
rehydrated in a humid chamber for 2 hours, baked in a 
dry 80° vacuum oven for 2 hours, rinsed to remove un- 
15 absorbed DNA and then treated with succinic anhydride 
to reduce non-specific adsorption of the labeled 
hybridization probe to the poly-l-lysine coated glass 
surface. Immediately prior to use, the immobilized DNA 
on the array was denatured in distilled water at 90° 
20 for 2 minutes. 

For the pooled chromosome experiment, the 16 
chromosomes of Saccharomyces cerevisiae were separated 
in a CHEF agarose gel apparatus (Biorad, Richmond, CA) . 
The six largest chromosomes were isolated in one gel 
25 slice and the smallest 10 chromosomes in a second gel 
slice. The DNA was recovered using a gel extraction 
kit (Qiagen, Chatsworth, CA) . The two chromosome pools 
were randomly amplified in a manner similar to that 
used for the target lambda clones. Following 
30 amplification, 5 micrograms of each of the amplified 

chromosome pools were separately random-primer labeled 
using Klenow polymerase (Amersham, Arlington Heights, 
IL) with a lissamine conjugated nucleotide analog 
(Dupont NEN, Bost n, MA) for the pool containing the 
35 six largest chromosomes, and with a fluorescein 
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conjugated nucleotide analog (BMB) for the pool 
containing smallest ten chromosomes. The two pools 
were mixed and concentrated using an ultrafiltration 
device (Am icon, Danvers, MA). 
5 Five micrograms of the hybridization probe 

consisting of both chromosome pools in 7.5 /il of TE was 
denatured in a boiling water bath and then snap cooled 
on ice. 2.5 /il of concentrated hybridization solution 
(5 x SSC and 0.1% SDS) was added and all 10 pi 

10 transferred to the array surface, covered with a cover 
slip, placed in a custom-built single-slide humidity 
chamber and incubated at 60° for 12 hours. The slides 
were then rinsed at room temperature in 0.1 x SSC and 
0.1%SDS for 5 minutes, cover slipped and scanned. 

15 A custom built laser fluorescent scanner was used 

to detect the two-color hybridization signals from the 
1.8 x 1.8 cm array at 20 micron resolution. The 
scanned image was gridded and analyzed using custom 
image analysis software. After correcting for optical 

20 crosstalk between the fluorophores due to their 
overlapping emission spectra, the red and green 
hybridization values for each clone on the array were 
correlated to the known physical map position of the 
clone resulting in a computer-generated color karyotype 

25 of the yeast genome. 

Figure 6 shows the hybridization pattern of the 
two chromosome pools. A red signal indicates that the 
lambda clone on the array surface contains a cloned 
genomic DNA segment from one of the largest six yeast 

30 chromosomes. A green signal indicates that the lambda 
clone insert comes from one of the smallest ten yeast 
chromosomes. Orange signals indicate repetitive 
sequences which cross hybridized to both chromosome 
pools. Control spots on the array confirm that the 

35 hybridization is specific and reproducible. 
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The physical nap locations of the genomic DNA 
fragments contained in each of the clones used as array 
elements have been previously determined by Olson and 
co-workers (Riles, et al.) allowing for the automatic 
5 generation of the color karyotype shown in Figure 7. 
The color of a chromosomal section on the karyotype 
corresponds to the color of the array element 
containing the clone from that section. The black 
regions of the karyotype represent false negative dark 

10 spots on the array (10%) or regions of the genome not 
covered by the Olson clone library (90%) . Note that 
the largest six chromosomes are mainly red while the 
smallest ten chromosomes are mainly green matching the 
original CHEF gel isolation of the hybridization probe. 

15 Areas of the red chromosomes containing green spots and 
vice-versa are probably due to spurious sample tracking 
errors in the formation of the original library and in 
the amplification and spotting procedures. 

The yeast genome arrays have also been probed with 

20 individual clones or pools of clones that are 

fluorescently labeled for physical mapping purposes. 
The hybridization signals of these clones to the array 
were translated into a position on the physical map of 
yeast. 

25 

Example 2 

Total cDNA Hybridized to Micro Arrays of 
cDNA Clones with Two-Color 
Fluorescent Detection 

30 24 clones containing cDNA inserts from the plant 

Arabidopsis were amplified using PCR. Salt was added 
to the purified PCR products to a final concentration 
of 3 x SSC. The cDNA clones were spotted on poly-1- 
lysine coated microscope slides in a manner similar to 

35 Exampl 1. Among the cDNA clones was a clone 
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representing a transcription factor HAT 4, which had 
previously been used to create a transgenic line of the 
plant Arabidopsis, in which this gene is present at ten 
times the level found in wild-type Arabidopsis (Schena, 
5 et al., 1992). 

Total poly-A xdRNA from wild type Arabidopsis was 
isolated using standard methods (Maniatis , et al. , 
1989) and reverse transcribed into total cDNA, using 
fluorescein nucleotide analog to label the cDNA product 

10 (green fluorescence) . A similar procedure was 

performed with the transgenic line of Arabidopsis where 
the transcription factor HAT4 was inserted into the 
genome using standard gene transfer protocols. cDNA 
copies of mRNA from the transgenic plant are labeled 

15 with a lissamine nucleotide analog (red fluorescence) . 
Two micrograms of the cDNA products from each type of 
plant were pooled together and hybridized to the cDNA 
clone array in a 10 microliter hybridization reaction 
in a manner similar to Example 1. Rinsing and 

20 detection of hybridization was also performed in a 

manner similar to Example 1. Fig. 8 show the resulting 
hybridization pattern of the array. 

Genes equally expressed in wild type and the 
transgenic Arabidopsis appeared yellow due to equal 

25 contributions of the green and red fluorescence to the 
final signal. The dots are different intensities of 
yellow indicating various levels of gene expression. 
The cDNA clone representing the transcription factor 
HAT4, expressed in the transgenic line of Arabidopsis 

30 but not detectably expressed in wild type Arabidopsis, 
appears as a red dot (with the arrow pointing to it) , 
indicating the preferential expression of the 
transcription factor in the red-labeled transgenic 
Arabidopsis and the r lative lack of expression of the 



WO 95/35505 



PCI7US95/07659 



37 

transcription factor in the gr en-labeled wild type 
Arabidopsis . 

An advantage of the microarray hybridization 
format for gene expression studies is the high partial 
5 concentration of each cDNA species achievable in the 10 
microliter hybridization reaction. This high partial 
concentration allows for detection of rare transcripts 
without the need for PCR amplification of the 
hybridization probe which may bias the true genetic 

10 representation of each discrete cDNA species. 

Gene expression studies such as these can be used 
for genomics research to discover which genes are 
expressed in which cell types, disease states, 
development states or environmental conditions. Gene 

15 expression studies can also be used for diagnosis of 
disease by empirically correlating gene expression 
patterns to disease states. 

Example 3 

20 Multiplexed Colorimetric H ybridization on 

? <?r* dded solid Support 
A sheet of plastic-backed nitrocellulose was 
gridded with barrier elements made from silicone rubber 
according to the description in Section IV-A. The 

25 sheet was soaked in 10 x SSC and allowed to dry. As 

shown in Fig. 12, 192 M13 clones each with a different 
yeast inserts were arrayed 400 microns apart in four 
quadrants of the solid support using the automated 
device described in Section III. The bottom left 

30 quadrant served as a negative control for hybridization 
while each of the other three quadrants was hybridized 
simultaneously with a different oligonucleotide using 
the open-face hybridization technology described in 
Section IV-A. The first two and last four elements of 
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each array are positive controls for the coiorimetric 
detection step. 

The oligonucleotides were labeled with fluorescein 
which was detected using an anti-f luorescein antibody 
5 conjugated to alkaline phosphatase that precipitated an 
NBT/BCIP dye on the solid support (Amersham) . Perfect 
matches between the labeled oligos and the M13 clones 
resulted in dark spots visible to the naked eye and 
detected using an optical scanner (HP ScanJet II) 

10 attached to a personal computer. The hybridization 
patterns are different in every quadrant indicating 
that each oligo found several unique M13 clones from 
among the 192 with a perfect sequence match. Note that 
the open capillary printing tip leaves detectable 

15 dimples on the nitrocellulose which can be used to 
automatically align and analyze the images. 

Although the invention has been described with 
respect to specific embodiments and methods, it will be 
20 clear that various changes and modification may be made 
without departing from the invention. 
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IT IS CLAIMED: 



1. A method of forming a microarray of analyte- 
assay regions on a solid support, where each region in 
5 the array has a known amount of a selected, analyte- 
specific reagent, said method comprising, 

(a) loading a solution of a selected analyte- 
specif ic reagent in a reagent-dispensing device having 
an elongate capillary channel (i) formed by spaced- 

10 apart, coextensive elongate members, (ii) adapted to 
hold a quantity of the reagent solution and (iii) 
having a tip region at which aqueous solution in the 
channel forms a meniscus, 

(b) tapping the tip of the dispensing device 

15 against a solid support at a defined position on the 
surface, with an impulse effective to break the 
meniscus in the capillary channel and deposit a 
selected volume of solution on the surface, and 

(c) repeating steps (a) and (b) until said array 
20 is formed. 



2. The method of claim 1, wherein said tapping is 
carried out with an impulse effective to deposit a 
selected volume in the volume range between 0.01 to 100 

25 nl. 

3. The method of claim 1, wherein said channel is 
formed by a pair of spaced-apart tapered elements. 

30 4. The method of claim l,. for forming a plurality 

of such arrays, wherein step (b) is applied to a 
selected position on each of a plurality of solid 
supports at each repeat cycle proceeding step (c) . 
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5. The method of claim 1, which further includes, 
after performing steps (a) and (b) at least one time, 
reloading the reagent-dispensing device with a new 
reagent solution by the steps of (i) dipping the 
5 capillary channel of the device in a wash solution, 
(ii) removing wash solution drawn into the capillary 
channel, and (iii) dipping the capillary channel into 
the new reagent solution. 

10 6. Automated apparatus for forming a microarray 

of analyte-assay regions on a plurality of solid 
supports, where each region in the array has a known 
amount of a selected, analyte-specif ic reagent, said 
apparatus comprising 

15 (a) a holder for holding, at known positions, a 

plurality of planar supports, 

(b) a reagent dispensing device having ah open 
capillary channel (i) formed by spaced-apart , 
coextensive elongate members (ii) adapted to hold a 

20 quantity of the reagent solution and (iii) having a tip 
region at which aqueous solution in the channel forms a 
meniscus, 

(c) positioning means for positioning the 
dispensing device at a selected array position with 

25 respect to a support in said holder, 

(d) dispensing means for moving the device into 
tapping engagement against a support with a selected 
impulse, when the device is positioned at a defined 
array position with respect to that support, with an 

30 impulse effective to break the meniscus of liquid in 

the capillary channel and deposit a selected volume of 
solution on the surface, and 

(e) control means for controlling said positioning 
and dispensing means. 
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7. The apparatus of claim 6, vh rein said 
dispensing means is effective to move said dispensing 
device against a support with an impulse effective to 
deposit a selected volume in the volume range between 

5 0.01 to 100 nl. 

8. The apparatus of claim 6, wherein said channel 
is formed by a pair of spaced-apart tapered elements, 

10 9. The apparatus of claim 6, wherein the control 

means operates to (i) place the dispensing device at a 
loading station, (ii) move the capillary channel in the 
device into a selected reagent at the loading station, 
to load the dispensing device with the reagent, and 

15 (iii) dispense the reagent at a defined array position 
on each of the supports on said holder. 

10. The apparatus of claim 6, wherein the control 
device further operates, at the end of a dispensing 

20 cycle, to wash the dispensing device by (i) placing the 
dispensing device at a washing station, (ii) moving the 
capillary channel in the device into a wash fluid, to 
load the dispensing device with the fluid, and (iii) 
remove the wash fluid prior to loading the dispensing 

25 device with a fresh selected reagent. 

11. The apparatus of claim 6, wherein said device 
is one of a plurality of such devices which are carried 
on the arm for dispensing different analyte assay 

30 reagents at selected spaced array positions. 

12. A substrate with a surface having a 
microarray of at least 10 3 distinct polynucleotide or 
polypeptide biopolymers per 1 cm 2 surface area, each 
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distinct biopolymer sample (i) being disposed at a 
separate, defined position in said array, (ii) having a 
length of at least 50 subunits, and (iii) being present 
in a defined amount between about 0.1 femtomole and 100 
5 nanomoles. 

13. The substrate of claim 12, wherein said 
surface is glass slide coated with poly lysine, and said 
biopolymers are polynucleotides. 

10 

14. The substrate of claim 12, wherein said 
substrate has a water- impermeable backing, a water- 
permeable film formed on the backing, and a grid formed 
on the film, where said grid (i) is composed of 

15 intersecting water-impervious grid elements extending 

from said backing to positions raised above the surface 
of said film, and (ii) partitions the film into a 
plurality of water-impervious cells, where each cell 
contains such a biopolymer array. 

20 

15. A substrate with a surface array of sample- 
receiving cells, comprising 

a water- impermeable backing, 

a water-permeable film formed on the backing, and 
25 a grid formed on the film, said grid being composed of 
intersecting water-impervious grid elements extending 
from said backing to positions raised above the surface 
of said film. 

30 16 • The substrate of claim 15, wherein the cells 

of the array each contain an array of biopolymers. 
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17. A substrate for use in detecting binding of 
labeled biopolymers to one or more of a plurality 
distinct polynucleotides, comprising 
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a non-porous, glass substrate, 

a coating of a cationic polymer on said substrate, 

and 

an array of distinct polynucleotides to said 
5 coating, where each biopolymer is disposed at a 
separate, defined position in a surface array of 
biopolymers. 

18. A method of detecting differential expression 
10 of each of a plurality of genes in a first cell type 
with respect to expression of the same genes in a 
second cell types, said method comprising 

producing fluorescence- labeled cDNA's from mRNA's 
isolated from the two cells types, where the cDNA's 
15 from the first and second cells are labeled with first 
and second different fluorescent reporters, 

adding a mixture of the labeled cDNA's from the 
two cell types to an array of polynucleotides 
representing a plurality of known genes derived from 
20 the two cell types, under conditions that result in 

hybridization of the cDNA's to complementary-sequence 
polynucleotides in the array; and 

examining the array by fluorescence under 
fluorescence excitation conditions in which (i) 
25 polynucleotides in. the array that are hybridized 

predominantly to cDNA's derived from one of the first 
and second cell types give a distinct first or second 
fluorescence emission color, respectively, and (ii) 
polynucleotides in the array that are hybridized to 
30 substantially equal numbers of cDNA's derived from the 
first and second cell types give a distinct combined 
fluorescence emission color, respectively, 

wherein the relative expression of known genes in 
the two cell types can b determined by the observed 
35 fluorescence emission color of each spot. 
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19. The method of claim 18, wherein the array of 
polynucleotides is formed on a substrate with a surface 
having an array of at least 10 2 distinct polynucleotide 
or polypeptide biopolymers in a surface area of less 

5 than about 1 cm 2 , each distinct biopolymer (i) being 

disposed at a separate, defined position in said array, 
(ii) having a length of at least 50 subunits, and (iii) 
being present in a defined amount between about .1 
femtomole and 100 nmoles. 

0 

20. The method of claim 19, wherein said surface 
is a glass slide coated with poly lysine, and said 
biopolymers are polynucleotides non-covalently bound to 
said poly lysine. 
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Use of a cDNA microarray to 
analyse gene expression 
patterns in human cane r 

Joseph DeRisi u , Loiita Penland : & 
Patrick O. Brown 2 (Group I ); 
Michael L Bittner 3 *, Paul S. Meltzer\ 
Michael Rax- 3 , Vidong Chen\ Van A. Su 3 & 
Jeffrey M. Trent 5 (Group 2) 



The development and progression of cancer 1 " 3 and 
the experimental reversal of tumorigenicity 4 - 5 are 
accompanied by complex changes in patterns of 
gene expression. Microarrays of cDNA provide a 
powerful tool for studying these complex phenom- 
ena 6 " 8 . The tumorigenic properties of a human 
melanoma cell line, UACC-903. can be suppressed 
by introduction of a normal human chromosome 6 ( 
resulting in a reduction of growth rate, restoration 
of contact inhibition, and suppression of both soft 
agar clonogenicity and tumorigenicity in nude 
mice 4 - 5 - 9 . We used a high density microarray of 
1,161 DNA elements to search for differences in 
gene expression associated with tumour suppres- 
sion in this system. Fluorescent probes for 
hybridization were derived from two sources of cel- 
lular mRNA [UACC-903 and UACC-903(+6)] which 
were labelled with different fiuors to provide a direct 
and internally controlled comparison of the mRNA 
levels corresponding to each arrayed gene. The flu- 
orescence signals representing hybridization to 
each arrayed gene were analysed to determine the 
relative abundance in the two samples of mRNAs 
corresponding to each gene. Previously unrecog- 
nized alterations in the expression of specific genes 
provide leads for further investigation of the genet- 
ic basis of the tumorigenic phenotype of these cells. 

DNA microarrays. containing UM total elements, 
including N70 different cPNAsand controls 1 *- 1 1 (see 
Methods), were printed robotieally onto a glass micro- 
scope slide in four quadrants covering an area of about 
I enr (Fig. I ). We prepared fluorescent cDNA probes 
using total poly (A)* mRNA from UACC-Vt)3 cells and 
UACC-903(+6) cells by labelling with a u recti and red 
fluor, respectively. A mixture of the two flourescently 
labelled probes was hybridized to the DNA microarray. 
This comparative hybridization method, coupled with 
the doping of synthetic standards and an estimation of 
statistically significant deviation for local background 
variance allowed a direct and quantitative comparison 
of the relative abundance of individual DNA sequences 
in this complex sample*"* We added a set of synthetic 
poly (A) + -tailed 'mRNAs' to the purified mRNA from 
each cell line as internal standards to assist in quantita- 
tion and estimation f experimental variation intro- 
duced during labelling and reading. Targets 
complementary to these standards were included, in 
duplicate, on the micmarray. Based on these standards. 
mRNA species comprising 1:10.01X1 of the mass of the 
poly ( AP RNA could readily be delected. 
In a representative two-colour fluorescent scan of all 



spond t genes preferentially expressed in the tumori- 
genic UACC-903 cell line, and the reddish spots corre- 
spond to genes preferentially expressed in the 
non-tumorigenic UACC-903(+6) cell line. Genes 
expressed at approximately equal levels in the two cell 
lines appear yellow or brown. A portion of the array at 
higher magnification highlights the diverse pattern of 
differential expression observed (Fig. 2b). In Fig. 2c rec- 
tangles corresponding to specific array elements are 
coloured to reproduce the hue and intensity of the fluo- 
rescent signal at each element. The hybridization signals 
from a duplicated set of genes are shown juxtaposed, to 
illustrate the reproducibility of the hybridization signals 
for each gene. 

To address the possibility that an apparent difference 
in expression might result from experimental variables 
unrelated to the difference in chromosomal composi- 
tion between the two cell lines, we examined the vari- 
ance in expression for 90 'housekeeping' genes. We 
selected these genes based on the assumption that they 
would not be differentially expressed between the two 
cell lines. The averaged red/green ratio for this subset of 
genes was 1.13. The averaged red/green ratio for the set 
of five internal standards was 0.97 { // = 10). The vari- 
ability in the expression level of the housekeeping genes 
probably overestimates the experimental variability in 
measuring differential expression. As a conservative stan- 
dard, an absolute fluorescent signal (red or green) with 
an intensity greater than that observed at the control 
array elements containing total human genomic DNA 
was considered to represent specific hybridization. Gene- 
specific hybridization was therefore only considered sig- 
nificantly different between samples if the following two 
criteria were met: i) the signal intensity (green or red) 
exceeded this threshold; and ii) the logarithm of the 
red/green fluorescence signal ratio differed by >3 S.D. 
Irom the mean logarithm of this ratio for the 'house- 
keeping' gene panel (that is. ratios <0.52 or >2.4). 

By these criteria, mRNA levels for 15/870 (1.7%) genes 
were significantly diminished, while the mRNA levels 
lor 63/870 (7.3%) genes were significantly increased in 
association with suppression of tumorigenicity by intro- 
duction of chromosome 6. To test the reliability of 
microarray hybridization results in identifying differen- 
tially expressed genes, we analysed 16 genes by north- 
ern analysis. In each case, the results of northern analysis 
corroborated the differential gene expression identified 
by microarray hybridization (Fig. 3). 

Significant differences in expression between these 
two cell lines identified several genes as candidates for 
determining features ol the tumorigenic phenotype of 
the melanoma cells. For example, among the genes 
delected with significantly higher expression (> 10-fold) 
in the tumorigenic cells was the human brown locus pro- 
tein (TRP I /melanoma antigen gp75). This is the most 
abundant glycoprotein in melanocyte cells and a critical 
melanosome membrane protein 12 -' \ Additionally, its 
expression is reduced when melanoma cell lines are 
induced to differentiate by treatment with HMBA 12 * 13 . 
Also expressed at a significantly higher level was a spliced 
variant of the mRNA encoding myelin IM.IVn.M20. This 
is widely expressed in neural crest derived cells in early 
development and has been suggested to play a role in 
cell-cell signaling during development 14 . 
• • J r ' — u - '• • *• »- : r — — ...i^*A lev- 
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Fig. 1 Properties of cONA rmooarTays. a. A fluorescent scan of DNA printed onto a pory-lysme coated slide. The DNA is stained with a DNA-specrlic Ituorescent 
dye. YOYO. The center-to-center spacing of adjacent spots is 450 U. allowing the potential for up to 10.000 spots/2.54 X 7.62 cm microscope slide. 6. Effi- 
cient blocking of hybridization to DNA repeats. Hybridization of fluorescein- labelled poly (dT)* to arrays in the absence of competitor produces strong 
hybridization t 0 immobilized poly (dA)* as well as to some cONAs. such as the EST T64827 shown. Rhodamine-labeHed cDNA (red) from the UACC-903 ceU 
lme hybridized m the presence of poly (dA)' blocker shows Mile it any signal at either site (Total H = total human). Similarly, hybndaauon with tluorescetn-tabelled 
Con DNA in the absence ol competitor produces bright signal on immobilized Cot1 DNA. total human DNA and at some cDNA elements (presumed to con- 
tain highly repeated sequences, such as R23416): while Rhodamine- labelled cDNA (red) from the UACC-903 cell line produces little it any signal at these 
locations when hybridized m the presence of excess untabelled poly (dA)*. and human Coll DNA. The absence of signal at some cDNA locations following 
UACC-903 cDNA hybridizations arso indicates that the PCR-amplified. ptasmid vector sequences at all cDNA targets oo not contribute significant hybridiza- 
tion signal, c. Schematic of the array organisation. Robotic printing from 96 well microliter trays was carried out with 4 print heads, spaced to tit into 4 adja- 
cent microliter wells. This maps the contents of each tray into (our separate quadrants on the glass slide. A colour-coded map ot the general distribution of 
target types m each of the resulting quadrants is shown. 



els were elevated by the addition of" a normal chromo- 
some 6(17 genes) are known to be activated by iFN-y. a 
cardinal proinflammatory cytokine that, a moot; other 
activities, induces expression of the gene products of the 
MHC class II locus. For example, the mRNA encoding 
monocyte chemotactic protein I (MCAF/MCF1 ), a 
cytokine that induces monocyte chemotaxis and aciiva- 
lion ,r "\ was more than 10- fold less abundant in the 
tumorigenic cell line. In the skin, MPCl is critical in the 
regulation of cutaneous monocyte trafficking" 1 " 1 ', and 
elevated expression plays a role in suppression of tumour 
growth and metastasis' 9 " 21 . The mechanism by which 
these interferon-y regulated genes are induced in UACC- 
903 cells by transfer of a normal chromosome 6 remains 
to be determined. It is worth noting, however, that the 
interferon-y receptor gene is localized to the distal long 
arm of human chromosome 6. 

Finally, several genes that showed > 10-fold higher 
expression in the suppressed UACC-903(+6) cells have 
previously been recognized in other models of tumour 
suppression. Most notably, there was elevated expres- 
sion of the mRNA encoding WAFl (p21 ), a key media- 
t r f tumour suppressi n by p53 (ref. 18). The p21 
protein had previously been identified as a melan ma 
differential n-associated antigen (termed mda-6) l9J0 . 
In melanoma cell lines suppressed for metastasis by the 
introduction of chromos me 6, express! n of WAFl 
(p21 ) mRNA and protein correlates inversely with 
metastatic potential 20 . 



These results provide a wide view of the diverse sys- 
tems that are altered in this model system of tumori- 
gen icily, and focus attention on specific gene products 
and pathways that may be of particular importance in 
this tumour type. 

Our ability to classify human cancers in a way that 
reflects the underlying molecular pathology or that 
anticipates their potential for progression or response 
to treatment, remains primitive. Using cDNA microar- 
rays to define alterations in gene expression associated 
with a specific cancer may be an efficient way to uncov- 
er clues to the specific molecular derangements that con- 
tribute to its pathogenesis and thus identify potential 
targets for therapeutic intervention. Moreover, recogni- 
tion of pathognomonic alterations in gene expression 
might provide a basis for improved diagnosis and mol- 
ecular classification of cancers and thus allow selection of 
the most appropriate therapeutic strategies. 

Public databases of human expressed gene sequences 
contain partial sequences of at least 40,000 different 
human genes 11 , and efforts to develop a human tran- 
script map have developed rapidly 21 . Based on the high 
yield of information obtained using an array of < 1 ,000 
different genes, a more comprehensive survey of gene 
expression patterns, using a more complete array of 
human genes, will likely pr vide a rich source of new 
and useful insights into human biology and a deeper 
understanding of the gene pathways involved in the 
pathogenesis f cancer and other diseases. 
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Methods 

Generation of microarrays. hybridization, scanning The 
preparation ol coated microscope slides and subsequent robot- 
ic priming nl I >\'A was carried mil m a mannci similar m ihat 
described . Ilrieliy. prc-elcaned glass slides wen- treated with 
|M»ly-I.. lysine solution tStiim.il to form an adhesive surface lor 
printing. PCK products, purified In cth.mol purification, were 
reMispendcd in 3\ SSC. A ciiNtom built arraying robot picked 
up and deposited small volumes (-5 nanohters! of I )N'A onto 
the sluice Alter printing ihe slides were washed in a H2"<» SI >S 
solution. The remaining hound I )NA was denatured In sub- 
merging the slides in 95 "C distilled water lor 2 min followed In 
a brief wash with y.Vi„ ethanol. UNA was I'V cmsslinked in the 
slides (Sirjtjjsenc Siratalinker. Ml m|). To prevent non-specific 
probe binding, the slides were blocked In rins.ng in a solution* 
nl 70 niM succinic anhydride dissolved in 0. 1 M* hone acid pH 
K.o. containing 35% I -methyl- J-pvrrolidinone (Aldrich). 
Addiiiiin.il protocols and parts list pertaining to microarrav 
fabrication can he obtained from hitp://cmgm.sianford.edu/ 
phrown. 

l J urilied. labelled cDNA was resiispendcd in I I pi of 3.5x SSCI 
containing 4 ug o| poly (dAK UNA. 15 "g K ioh (UNA, 4 ^ol" 
human Cm I I )NA (Cibco ItKU. and 0.3 ul of" 10",. SI )S. Prior to 
hybridization, the .solution was boiled lor : nun then allowed to 
cool lo i. mm lein|H*rature. I lyhridi/alion was carried out ai 
tO °C lor -l-l h in a water bath. Prior to scanning, .slides woe 
washed in 2* Sm ;. ().:<■.., S!>S lor 5 min and 0..V S^( '. Un I min, 

. - — 

I.kmI l.iser m.er..>eope huili hv S. Smith with software wriiien 



by N. A separate scan, iimh- the appropriate excitation line, 
u-jx done lor eaeh oi the tu,» llutuophores used. Data was col- 
lected .ii .i maximum resolution „fy nnei*»n.s/pixel with 12 bits 
ol depth. 

Probe preparation and labelling. UNA was extracted from cells 
usin^ ilu'Tiia/ol reagent < I I I Inc.). followini: the maiuilactur- 
ei s directions, c I )N A pnibes were synthesi/ed from singly oltgo 
il l-selected (Pharmac.al mKNA pool.s. Huorescentlv bMled 
cDNA was prepared from mU\A by olij;*. d*| - primed |H)lvmer- 
i/ai.on usin- Superscript II reverse t'ranscripuse (LTl 'lnc). 
The pool of nucleotides in the labelling reaction was 0.5 mM 
dCTP. dATP and dC I P and ().: mM JT|'|» Fluorescent 
nucleotides, khodaminc III) dUTP fPerkin F.lmer Ctflus) or 
Cv3 dlTP I Amer.shanij. were preNent at 0.1 mM. Probes were 
purified hv gel chromatography (bioSpin 6/HioUad) and 
elhanol precipitation. 

Selection of cDNA elements and generation of control tern- 
plates. Synthetic cDNAs were prepared by cloning random 
/in/i/HI and Hmdlll ended fragments of E. toli DNA in the vec- 
tor pSPM poly (A)' (Promega). linearizing iMilated plasmid 
DNA with K(f»lil and synihesi/ing p.»|y (A)* tailed RNA com- 
plementary to the insert horn the resident SPn promoter 
< Piomega). Prior to use. the >\ nihesi/ed liNAs were xHixied on 
oligo d I cellulose.' The largest group ol cDNAs consisted of 674 
cDNA climes from the INIh arraved normali/ed tnt'ant brain 

• * M.^ieJ u» i; ^. .-.ere MB 

library member thai on responded to a named gene according 
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MARCKS 



Fig. 3 Nonnern nyoridizaien substantiat- 
ing the consistency I the cDNA microar- 
ray resutts. Corresponding locations within 
the cDNA microarray illustrated m Fig. 2a 
are provided (or 1) Waf-Vp2V. 2) MARCKS: 
3) collagenase: 4) MCAF/MCP-1; 5) n- J- 
anttchymotrypsin: and 6) 0 -actm. The sig- 
nal detected by a radio-labelled p-actm 
probe represents a control for loading vari- 
ance, with a red/green ratio observed on 
the cDNA microarray (fig. 2a.c) for f^actin 
of V04. 




Collagenase 



MCAF/MCP-1 



to the UniCcnc EST clustering sys- 
lcm :, - : . The second largest group of 
clones consisted of" 1 83 sequenced cDNA 
cloneN genera led by subtraction of cDNA 
from the chronioMime-o suppressed 
non-tumorigenic L'ACC-9(»3 (+6) cell 
line with c I >NA from its parental htinort- 
genie cell line UACC-903 (ret. V). 
Approximately KM) additional genes 
(total H70 genes arrayed) were obtained 
from KST libraries on the basis of their 
r> 7 -Antichymotrypsin ex press i< »n pal tern 1 1 issue specific, and so 
on), l-uh array included the following 
hybridization controls: plasm id vector, 
lambda. OX 17-1 phage, total human DNA. 
human Gill l>NA. and poly (A)\ The 
synthetic standards used for normal iza- 
tion of signals in each wavelength were 
also arrayed. Gmtrols were included in 
each ijtiadrani of the array lo assess the reprodticibility of the 
hybridization signal. Two plates of cDNA clones (derived from 
the UAGi-W subtracted library) were also arrayed in dupli- 
cate. Tideliiy of the Unigene array relative to dhKST was tesied 
by sequencing of a random sample of 1 1 clones used lor 
microarray c* in st ruction. All .sequences wvre identical with the 



P-Actin 



corresponding dhl.ST entries. Additionally, each muinjii.i U j 
cDNA from the LWlX-Wtt Mibtracioi library wji H-queiued. \ 
listing ol cDNAs comprising this miaoarr.ty %%huh woe 
derived I rum the Unigene and 'housekeeping* panel fc an Iv 
obtained from hltp://w\w.nih.gi»v/|)IU/LCr./ARKAVc\pn - 
html. 

Northern blot analysis Total RNA, 10 ug per lane, was elec* 
irophoresed in 1.2% agarose- formaldehyde gels and transferred 
onto nylon membrane (Hybond-N*. Amcrsham) by capitlary 
blotting overnight. For DNA probes insert fragments from the 
Snares 1Mb cDNA library'" were obtained by vector PCR for 
p2l. MAKCKS, U- l-amichymotrypsin and (5-actin. Probes for 
fibroblast collagenase and MCAF/Ml^rM were isolated from a 
UACC-903t+fc) enriched cDNA library** with all probes 
labelled by random priming. Fillers were washed to a mi in- 
gency oft). I x NSC at -12 *C for 20 min. 

Web sites, http://cmgm.stanlord.edu/pbrown for protocols and 
parts list pertaining to microarray fabrication. 
http://www.nchgr.nih.gOv/l )IR/LCi;/AURAY/expn.html for a 
listing til cPNAs comprising ibis microarray which were 
derived from the Unigene and 'housekeeping panel. 
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A DNA Microarray System for Analyzing " 
Complex DNA Samples Using Two-color 
Fluorescent Probe Hybridization 

Dari Shalon, 1 - 4 Stephen J. Smith, 3 and Patrick O. Brown 1 ' 2 ' 5 

toward Hughes Medical Institute and Departments of 2 Biochemistry and J Molecular and Cellular 
Physiology, Stanford University, Stanford, California 94305 



Detecting and determining the relative abundance of diverse individual sequences in complex DNA samples is 
a recurring experimental challenge in analyzing genomes. We describe a general experimental approach to 
this problem, using microscopic arrays of DNA fragments on glass substrates for differential hybridization 
analysis of fluorescently labeled DNA samples. To test the system, 864 physically mapped X clones of yeast 
genomic DNA, together representing >7S% of the yeast genome, were arranged into 1.8-cm x 1.8<m arrays, 
each containing a total of 1744 elements. The microarrays were characterized by simultaneous hybridization 
of two different sets of isolated yeast chromosomes labeled with two different fluorophores. A laser 
fluorescent scanner was used to detect the hybridization signals from the two fluorophores. The results 
demonstrate the utility of DNA microarrays in the analysis of complex DNA samples. This system should 
find numerous applications in genome-wide genetic mapping, physical mapping, and gene expression studies. 



Many problems in genome analysis depend on 
determining what specific sequences are repre- 
sented in a complex DNA or RNA sample and at 
what abundance, for example, what genes are 
represented in a specific chromosome band or 
YAC clone, what intervals are amplified or de- 
leted in a particular cancer cell, or what genes are 
expressed in specific cells under specific condi- 
tions. As a general approach to this problem, we 
have developed a system for making microarrays 
of DNA samples on glass substrates, probing 
them by hybridization with complex fluorescent- 
labeled probes, and using a laser-scanning micro- 
scope to detect the fluorescent signals represent- 
ing hybridization. Fluorescent labeling allows for 
simultaneous hybridization and separate detec- 
tion of the hybridization signal from two or more 
probes. This in turn allows very accurate and re- 
liable measurement of the relative abundance of 
specific sequences in two complex samples. 

RESULTS 

Array Hybridization Pattern 

Figure 1 shows the two-color fluorescent scan of 
a yeast genomic array following hybridization 

'Preterit address: Syntenl, Inc., Palo Alto. California 94305. 
'Corresponding author. 

E-MAIl pbrown cm9m.stanf0rd.edu, http://cmgm. 
stanford.edu/pbrown; FAX (415) 72 J 1399. 



with a mixed probe consisting of lissamine- 
labeled DNA from the 6 largest yeast chromo- 
somes together with fluorescein-labeled DNA 
from the 10 smallest yeast chromosomes. A red 
color indicates that yeast sequences present in 
the lissamine-labeled hybridization probe hy- 
bridized to an array element. A yellow-green 
color indicates that yeast sequences present in 
the fluorescein-labeled hybridization probe hy- 
m bridized to an array element. An orange color in- 
dicates cross-hybridization of both chromosome 
pools to an. array element (e.g., dispersed repeti- 
tive elements, such as Tyl elements). 

Each clone was spotted twice, resulting in du- 
plicate hybridization patterns in adjacent quad- 
rants of the array. Control DNA spots, which 
were randomly amplified in the same manner as 
the X clone array elements, are located in the bot- 
tom corner of each quadrant. "A" points to a pair 
of spots containing total yeast genomic DNA. 
These spots appear orange because both chromo- 
some pools hybridized to yeast genomic DNA. 
The negative controls are as follows: "B" points 
to a pair of spots of wild-type X DNA, "C" points 
to a pair of human genomic DNA spots, and "D" 
points to a pair of 6X174 DNA spots. The lack of 
a hybridization signal at these three negative 
control spots indicates that the hybridization was 
specific for yeast sequences. 
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Figure 1 Two<olor fluorescent scan of a 1 .8-cm x 1 .8-cm yeast array 
of K clones of yeast genomic DNA. The DNA spots are spaced at a 
distance of 380 p.m from center to center. A probe mixture consisting of 
DNA f r m the 6 largest yeast chromosomes (4, 7, 1 2, 1 3, 1 5, 1 6) labeled 
with lissamine (red dots) and DNA from the 10 smallest yeast chromo- 
somes (1, 2, 3, 5, 6, 8, 9, 10, 11, 14) labeled with fluorescein (yellow- 
green dots) was hybridized to the array. A pair of yeast genomic DNA 
spots (A) served as a positive control. The three negative controls are X 
DNA (B), human genomic DNA (C), and 6X174 DNA (D). 



Karyorype Depiction of the Array Hybridization 
Pattern 

The inserts contained in the arrayed X clones 
have been mapped physically (Riles et al. 1993). 
The clones are arrayed in a random but known 
order on the array. Therefore, using the identity 
of each clone along with its physical map infor- 
mation, the pattern of hybridization to the yeast 
array can be represented in the form of a karyo- 
type of the yeast genome, as shown in Figure 2. 
The color of any segment of the ideogram repre- 
senting an individual chromosome on the karyo- 
type is directly determined by the ratio of red and 
green hybridization signals at the array positions 
of the corresponding clones. The lengths of the 
discrete colored segments of each chromosome 
correspond to the physical lengths of the yeast 



inserts. The chromosome seg- 
ments colored black represent ei- 
ther intervals of the genome that 
are not represented by clones in 
the library (90%) or false-negative 
hybridization signals on the array 
(10%i. Most of these false nega- 
tives are attributable to failures of 
the PCR amplification of the X 
clones, though occasional failures 
of the arraying process or nonuni- 
form surface preparation could ac- 
count for a small fraction of the 
false-negative signals. The large 
gap on chromosome 12 is the re- 
gion coding for ribosomal DNA 
that was not represented among 
the arrayed clones. Genomic inter- 
vals represented by overlapping 
clones were assigned a color based 
on the hybridization signals of 
only one of the overlapping 
clones, chosen at random. 

Note that in this representa- 
tion of a yeast karyotype, the larg- 
est six chromosomes are mainly 
colored red. This indicates that 
most of the arrayed clones that 
were mapped previously to these 
six large chromosomes hybridized 
primarily to the lissamine-labeled 
probe prepared from the corre- 
sponding purified chromosomes. 
Conversely, the smallest 10 chro- 
mosomes are mainly colored green 
in this image, matching the origi- 
nal CHEF gel isolation of the chro- 
mosomes used as the hybridization probe. The 
experiment was repeated with the yeast genome 
split into six discrete chromosome pools contain- 
ing 2-4 chromosomes per pool using CHEF gel 
electrophoresis. The chromosomes in each pool 
were extracted from the gel, amplified, and fluo- 
rescently labeled. The six chromosome pools 
were hybridized to six separate yeast arrays. 
Forty-four X clones gave a positive hybridization 
signal on all six arrays indicating that they con- 
tain yeast repetitive sequences (data not shown). 
These 44 clones and 10 clones with very weak 
hybridization signals were not included in the 
data set used to produce this karyotype. 

There were -40 anomalous clones, which ap- 
pear in this karyotype representation as green 
bands on the otherwise red chromosomes or red 
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Figure 2 Computer-generated ideogram repre- 
senting a karyotype of 5. cerevisiae, based on the 
normalized hybridization signals from the array 
shown in Fig. 1 . Note that the 6 largest chromo- 
somes are mainly red and the 10 smallest chromo- 
somes are mainly green. Black stripes represent in- 
tervals not represented by clones in the array or for 
which the corresponding clones gave false-negative 
hybridization signals. 

bands on the otherwise green chromosomes. 
Four randomly chosen examples of these anoma- 
lous clones were analyzed by hybridizing the 
clones to vertical strips cut from a Southern blot 
of CHEF gel-separated yeast chromosomes. In 
each case, the hybridization patterns of the 
anomalous clones corroborated the chromo- 
somal locations assigned by the microarray hy- 
bridization results (data not shown). Two clones 
that were thought to map to the 10 smallest chro- 
mosomes were found to hybridize preferentially 
to the probe representing the 6 largest chromo- 
somes and thus appear as anomalous red bands 
on the karyotype. Both hybridized to one of the 
six largest chromosomes on the Southern blot. 
Similarly, two clones that appear as anomalous 
green bands on the karyotype were found to hy- 
bridize to one of the 10 smallest chromosomes on 
the Southern blot. Thus, the anomalous clones 
are probably the result of sample tracking errors 
or, possibly, of errors in the published restriction- 
digest-based physical map on which the karyo- 
type representation was based (Riles et al. 1993). 

DISCUSSION 

The DNA microarray hybridization system re- 
ported here is conceptually and functionally 



similar to fluorescent in situ hybridization (FISH) 
to metaphase chromosomes, with three impor- 
tant differences. First, the target elements of the 
microarrays can, in principle, be any length or 
composition, from megabase YAC clones or mi- 
crodissected chromosome bands to individual 
cDNA clones, to short oligonucleotides. This ver- 
satility allows the user to choose characteristics, 
such as the mapping resolution and genetic com- 
plexity of each array element, to suit a particular 
application. Second, the hybridization signals are 
localized to discrete elements of known size and 
location, making them easier to identify and 
quantitate than the hybridization signals from 
irregularly shaped metaphase spreads. Third, mi- 
croarrays are more consistent and potentially 
amenable to automated production, hybridiza- 
tion, and data analysis than metaphase spreads. 

Arrays of DNA samples on porous mem- 
branes, for example, dot blots, have long been 
used as a basic tool in molecular biology. Dot- 
blot membranes are usually at least 8 x 12 cm in 
size, require the use of milliliter volumes of hy- 
bridization solution, and are limited, owing to 
autofluorescence and scattering, to radioactive, 
chemiluminescent, and colorimetric hybridiza- 
tion detection methods (Ross et al. 1992). Micro- 
arrays made on glass surfaces, on the other hand, 
can be mass-produced and are comparatively in- 
expensive, convenient, and compatible with 
fluorescent hybridization detection methods. 
Furthermore, a glass surface, when appropriately 
treated, has very low nonspecific binding of la- 
beled hybridization probes, resulting in lower 
backgrounds than are encountered typically with 
porous membranes. For hybridizations with very 
complex probes, the concentration of the labeled 
probe DNA is a limiting factor in the sensitivity 
of the assay. Minimizing the volume of the probe 
solution in a hybridization, by restricting the tar- 
get to a small area and by using a nonporous 
substrate, makes it practical to achieve very high 
probe concentrations. 

One important advantage of fluorescently la- 
beled probes is that, unlike most radioactive and 
chemiluminescent signals, fluorescent signals do 
not disperse and therefore allow for very dense 
array spacing. A unique, and probably the m st 
important, advantage of fluorescent probes is 
that the hybridization signals from two or more 
differently labeled probes hybridized to the same 
target element can be detected separately. In this 
way, two-color hybridization detection allows for 
a direct and quantitative comparison f the 
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abundance of specific sequences between two 
probe mixtures that are hybridized competitively 
to a single array. The absolute intensity of a hy. 
bridization signal at a particular element in an 
array can vary owing to experimental factors 
such as variations in the amount of DNA depos- 
ited on the array, variations in the hybridization 
or wash conditions between experiments, or 
variations in the hybridization characteristics of 
the different DNA sequences on the array. The 
ratio of the two signals at any element in an ar- 
ray, however, is relatively insensitive to these 
confounding factors because they affect both 
probe mixtures equivalently. This ratio therefore 
accurately reflects the relative abundance of the 
cognate sequence in the two probe samples. This 
is the principle underlying the technique of com- 
parative genomic hybridization (CGH), which is 
used to detect changes in the copy number of 
specific chromosomes or chromosomal regions 
(Kallioniemi et al. 1992). CGH is based on mea- 
suring the relative fluorescent hybridization in- 
tensities of two genomic-complexity hybridiza- 
tion probes, for example, probes representing ge- 
nomic DNA from normal and affected tissue 
samples, which are labeled with two distinct fluo- 
rophores and hybridized simultaneously to a 
metaphase spread. DNA microarxay representa- 
tions of the human genome may provide a more 
convenient and higher resolution alternative to 
metaphase chromosomes for CGH. 

Cross-hybridization between related se- 
quences is an important problem faced by any 
hybridization-based assay, including the DNA 
microarray assay described here. Studies are now 
in progress to quantitate the extent of cross- 
hybridization between related sequences of vary- 
ing homology and length, in DNA microarray 
hybridizations. The stringency of hybridization 
and washing can be controlled by varying the salt 
concentration and temperature as in conven- 
tional membrane-based hybridizations. Cross- 
hybridization caused by repetitive sequences can 
be minimized by prehybridization of the probe or 
array with vast excess of unlabeled copies of the 
repetitive sequences. 

Alternative methods have been described for 
making microarrays of very short DNA se- 
quences, involving photolithography (Pease et 
al. 1994) or physical masking (Maskos and South- 
ern 1992) methods. These in situ synthesis meth- 
ods are inherently limited to low complexity ar- 
ray elements consisting of oligonucleotides. For 
complex-probe hybridizations, the specificity of 



hybridization is improved by using DNA frag- 
ments substantially longer than oligonucleo- 
tides. Moreover, the in situ synthesis approaches 
to array fabrication depend on pri r knowledge 
of the sequence to be recognized by each array 
element. The approach described here makes mi- 
croarrays by transferring tiny volumes of DNA 
samples from microwell storage plates to a solid 
substrate. Thus, nucleic acids (or other mol- 
ecules) of virtually any length or any origin can 
be arrayed, and knowledge of their sequences is 
not required. 

The arrays used in these experiments do not 
represent the maximal achievable density of ele- 
ments. We have found that the spacing between 
the spots can be decreased by shrinking the con- 
tact area of the printing tip and by increasing the 
hydrophobicity of the glass surface. Microarrays 
with 100-M.m feature size have been tested suc- 
cessfully in pilot experiments (data not shown). 
Assuming the projected availability of the appro- 
priate physically mapped human genomic clones 
(Hudson et al. 199S), arrays at 100-n.m spacing 
would allow for 10,000 discrete intervals of the 
human genome to be represented in a 1-cm 2 ar- 
ray. Such an array could be used for mapping at a 
resolution of <0.5 Mb. Experiments are in 
progress to explore the feasibility of such arrays. 

Our initial motivation for developing these 
microarrays arose from the need for abundant 
and inexpensive genomic arrays for genomic 
mismatch scanning (GMS) (Nelson et al. 1993), a 
method of genetic linkage analysis based on 
identification of the regions of "identity by de- 
scent" between affected relative pairs using a 
single complex-probe hybridization to an array 
of genomic clones. Experiments using these ar- 
rays to map quantitative trait loci in yeast by 
GMS are currently in progress Q. deRisi, D. Lash- 
kari, L Penland, L McAllister, J. McCusker, R. 
Davis, and P.O. Brown, unpubl.). 

Microarrays of cDNA clones, prepared using 
the system described here, have been used for 
quantitative monitoring of gene expression pat- 
terns in Arabidopsis (Schena et al. 1995), S. cercvi- 
siae (D. Lashkari, J. deRisi, L Penland, P.O. 
Brown, and R. Davis, unpubl.), and human tis- 
sues 0. deRisi, M. Bittner, P. Meltzer, L. Penland, 
J. Trent, and and P.O. Brown, unpubl.). We an- 
ticipate that DNA microarrays of the kind de- 
scribed here will be useful in additional applica- 
tions for which conventional dot blots, high- 
density gridded arrays on porous membranes, or 
FISH are currently used. These potential applica- 



642 ^ GENOME RESEARCH 



DNA MICRO ARRAYS FOR ANALYZING COMPLEX D\A SAMPLES 



tions include comparative genomic hybridiza- 
tion (KaUioniemi et al. 1992), sequencing by hy- 
bridization (Drmanac et al. 1993), physical map- 
ping of cloned or amplified sequences (Billings et 
al. 1991), and economical distribution of re- 
agents for integrated genetic and physical map- 
ping based on a common set of arrayed clones 
(Zehetner and Lehrach 1994). 

METHODS 

Amplification of Target DNA Elements 

The array elements were prepared from physically mapped 
x cl nes (Riles et al. 1993). The X clones were amplified 
using randomly primed polymerase chain reaction (PCR) 
based on published and unpublished protocols (Bohlander 
et al. 1992; S. Nelson, unpubl.). The phage lysates were 
amplified in a 10-*il PCR reaction using 5 »m final concen- 
tration of primer A (GCTATCTTCAAGATCANNNNNN), 
200 *lm dNTPs. and 1 unit of Taq polymerase. Round A 
consisted of five cycles at 94 # C for 1 min, 25*C for 1.5 min, 
25-72°C over 7 min, and 72°C for 3 min using Taq poly- 
merase (BMB). For round B, the reaction volume was 
brought up to ] 00 mJ for a final concentration of 2 of 
primer B (GCTATCTTCAAGATCA), 200 hm dNTPs, and 4 
units f Taq polymerase. Round B consisted of 30 cycles of 
94°C f r 1 min/ 56 B C for 2 min, and 72*C for 3 min. The 
amplification was performed in 96-well plates using crude 
phage lysates as the templates, resulting in an amplifica- 
tion of both the 35-kb X vector and the 5-kb to 15-kb yeast 
insert sequences as a distribution of PCR products between 
250 bp and 1500 bp in length. 

The PCR products were purified and transferred into 
TE (10 mM Tris, 1 mM EDTA at pH 8.0i buffer using Sepha- 
dex G50 gel filtration (Pharmacia) and evaporated to dry. 
ness at room temperature overnight. Each of the 864 am- 
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Figure 3 The layout of the arraying machine. All motions are under computer 
control. For more details of the arraying machine, see web page http:// 
cmgm.stanford.edu/pbrown. 



plified X clones was rehydrated in 15 nl of 3 x SSC 
(20 x SSC = 3 m NaG, 0.3 m Na 3 citrate) in preparation for 
spotting onto the glass under norma) room temperature 
conditions. 



Preparation of DNA Microarrays 

The microarrays were fabricated on poly-i-lysine coated 
microscope slides (Sigma). A custom-built arraying ma- 
chine, consisting of four tweezer-like printing tips 
mounted 9 mm apart on a computer-controlled robotic 
stage (Shalon 1996), loaded 1 p.1 of the concentrated PCR 
product directly from corresponding clusters of four wells 
of 96-well storage plates and deposited -5 nl of each 
sample onto each of 40 slides. Surface tension loaded the 
sample into the printing tip directly from the microwell 
plate and held the sample in the tip during the printing 
operation. Printing was achieved by lightly tapping the tip 
against the glass surface. The open-capillary design al- 
lowed for rapid rinsing and drying of the tips between 
samples. Figure 3 shows the layout of the arraying ma- 
chine. Figure 4 shows a detailed view of the four printing 
tips and the staggered printing partem on the microscope 
slides. Adjacent samples were spotted 380 n-m apart on the 
slides. After each set of four samples was printed onto 40 
slides, the printing tips were rinsed with a jet of water for 
2 sec and then dried by lowering the tips onto a sponge for 
2 sec. The process was repeated for all 864 samples and 
eight control spots. 

After the spotting operation was complete, the slides 
were rehydrated in a humid chamber at room temperature 
for 2 hr, baked in an 80*C vacuum oven for 2 hr, then 
rinsed in 0.1% sodium dodecyl sulfate (SDS) to remove 
un adsorbed DNA. To reduce nonspecific adsorption of the 
labeled hybridization probe to the poly-L-lysine coated 
glass surface, the slides were treated with succinic anhy- 
dride. One gram of succinic anhydride was dissolved in 
100 ml of l-methyl-2-pyrrolidinone and then 100 ml of 
0.2 m boric acid (pH 8.0) was 
added. The arrays were soaked in 
this solution for 10 min and then 
rinsed in distilled water four 
times for 5 min each. Immedi- 
ately before use, the arrayed DNA 
elements were denatured by plac- 
ing the slide in distilled water at 
90*C for 2 min. 



Amplification and Labeling 
of Hybridization Probe 

The 16 chromosomes of Saccharo- 
myees eerevisiae were separated us- 
ing a con tour -clamped homoge- 
neous electric field (CHEF) aga- 
rose gel apparatus (Bio-Rad) (Chu 
et al. 1986). The 6 largest chromo- 
somes were isolated in one gel 
slice and the smallest ten chro- 
mosomes in a second gel slice. 
The DNA from each slice was re- 
covered using a gel extraction kit 
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Figure 4 A close-up view of the four open- 
capillary printing tips. The tips are 9 mm apart and 
fit into four adjacent wells of a standard microwell 
plate and print arrays in a staggered fashion on mi- 
croscope slides. For more details of the printing tips, 
see web page http://cmgm.stanford.edu/pbrown. 



(Qiagen) and randomly amplified in a manner similar to 
that used in amplifying the target A clones (Grothues et al. 
1993). The main difference between this amplification 
procedure and the one used for the k array elements is a 
filtration step between rounds A and 8 to remove primer- 
dimers and the use of a random 9-mer 3' end on primer A. 
Following amplification, 2.5 jig of each of the amplified 
chromosome pools were separately random-primer labeled 
using Klenow polymerase (Amersham) with a lissamine- 
conjugated nucleotide analog (DuPont NEN) for the pool 
containing the 6 largest chromosomes and with a fluores- 
cein -conjugated nucleotide analog (BMB) for the pool con* 
taining the smallest 10 chromosomes. The two fluores- 
cent-labeled pools were mixed and concentrated using an 
ultrafiltration device (Amicon). 



Hybridization 

Five micrograms of the hybridization probe, consisting of 
both chromosome pools in 7.5 mJ of TE, was denatured in 
a boiling water bath and then snap-cooled on ice. Concen- 
trated hybridization solution (2.5 p.1) was added to a final 
concentration of 5x SSC/0.1% SDS. The entire 10 mJ of 
probe solution was transferred to the arcay surface, covered 
with a coversiip, placed in a custom-buiit single-slide hu- 
midity chamber, and incubated in a 60*C water bath for 12 
hr. The custom-built waterproof slide chamber has a cavity 
just slightly bigger than a microscope slide and was kept at 
100% humidity internally by the addition of 2 pj of water 
in a corner of the chamber. The slide was rinsed in 5 x 
SSC/0.1% SDS for 5 min and then in 0.2 x SSC/0.1% SDS 
for 5 min. All rinses were at room temperature. The array 
was then air dried, and a dr p of antifade (M lecular 
Probes) was applied to the array under a 24-mm x 30-mm 
coversiip in preparation f r scanning. 



Detection and Analysis 

A custom-built laser scanner was used to detect the two- 



color fluorescence hybridization signals from 1.8- 
cm x 1.8-cm arrays at 20-mjti resolution. The glass sub* 
strate slide was m unted n a computer-con trolled, two- 
axis translation stage (PM-500. Newport. Irvine. CA) that 
scanned the array over an upward-facing microscope ob- 
jective (20 x, 0.75NA Fluor, Nikon. Melville. NY) in a bi- 
directional raster pattern. A water<ooled Argon/Krypton 
laser (Innova 70 Spectrum. Coherent. Palo .Alto. CA). op- 
erated in multiline mode, allowed for simultaneous speci- 
men illumination at 488.0 nm and 568.2 nm. These two 
lines were isolated by a 488/568 dual-band excitation filter 
(Chroma Technology, Brattleboro. VT). An epifluores- 
cence configuration with a dual-band 488/568 primary 
beam splitter (Chroma) excited both fluorophores simul- 
taneously and directed fluorescence emissions toward the 
two-channel detector. Emissions were split by a secondary 
dichroic mirror with a 565 transition wavelength onto two 
multialkali cathode photomultiplier tubes (PMT; R928, 
Hamamatsu, Bridgewater. NJ), one with an HQ535/50 
bandpass barrier filter and the other with a D630/60 band- 
pass barrier filter (Chroma). Preamplified PMT signals were 
read into a personal computer using a 12-bit analog-to- 
digital conversion board (RTI-834, Analog Devices, Nor- 
wood, MA), displayed in a graphics window, and stored to 
disk for further rendering and analysis. The back aperture 
of the 20 x objective was deliberately underfilled by the 
illuminating laser beam to produce a large-diameter illu- 
minating spot at the specimen (5-p.m to 10-nm half- 
width). Stage scanning velocity was 100 mm/sec. and PMT 
signals were digitized at 100 *sec intervals. Two successive 
readings were summed for each pixel such that pixel spac- 
ing in the final image was 20 jun. Beam power at the 
specimen was -5 mW for each of the two lines. 

The scanned image was despeckled using a graphics 
program (Hijaak Graphics Suite) and then analyzed using 
a custom image gridding program that created a spread- 
sheet of the average red and green hybridization intensi- 
ties for each spot. The red and green hybridization inten- 
sities were corrected for optical cross talk between the fluo- 
rescein and lissamine channels, using experimentally 
determined coefficients. 
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